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Abstract 

We use the AdS/CFT correspondence to argue that large rotating black holes 
in global AdS^) spaces are dual to stationary solutions of the relativistic Navier- 
Stokes equations on S D ~ 2 . Reading off the equation of state of this fluid from the 
thermodynamics of non-rotating black holes, we proceed to construct the nonlinear 
spinning solutions of fluid mechanics that are dual to rotating black holes. In all 
known examples, the thermodynamics and the local stress tensor of our solutions 
are in precise agreement with the thermodynamics and boundary stress tensor of the 
spinning black holes. Our fluid dynamical description applies to large non-extremal 
black holes as well as a class of large non-supersymmetric extremal black holes, but 
is never valid for supersymmetric black holes. Our results yield predictions for the 
thermodynamics of all large black holes in all theories of gravity on AdS spaces, for 
example, string theory on AdS$ x S 5 and M theory on AdS 4, x S 7 and AdSj x S 4 . 
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1 Introduction 

In this paper, we predict certain universal features in the thermodynamics and other clas- 
sical properties of large rotating black holes in global AdSi) spaces for arbitrary D. Our 
analysis applies to black holes in any consistent theory of gravity that admits an AdS^ 
background; for example, IIB theory on AdS$ x S 5 or M theory on AdSj x S 4 or AdS± x S 7 . 

All theories of gravity on an AdS^ background are expected to admit a dual description 
as a quantum field theory on S D ~ 2 x time [H, 0]. Moreover, it is expected to be generally 
true that quantum field theories at sufficiently high energy density admit an effective 
description in terms of fluid dynamics. Putting together these facts, we propose that large, 
rotating black holes in arbitrary global AdSu spaces admit an accurate dual description 
as rotating, stationary configurations of a conformal fluid on S D ~ 2 . 

Assuming our proposal is indeed true, we are able to derive several properties of large 
rotating AdS black holes as follows: We first read off the thermodynamic equation of 
state of the dual 'fluid' from the properties of large, static, non-rotating AdS black holes. 
Inputting these equations of state into the Navier-Stokes equations, we are then able to 
deduce the thermodynamics of rotating black holes. In the rest of this introduction, we 
will describe our proposal and its consequences, including the tests it successfully passes, 
in more detail. 

Consider a theory of gravity coupled to a gauge field (based on a gauge group of 
rank c) on AdS/). In an appropriate limit, the boundary theory is effectively described 
by conformal fluid dynamics with c simultaneously conserved, mutually commuting U(l) 
charges Ri (i = 1 . . . c) . Conformal invariance and extensivity force the grand canonical 
partition function of this fluid to take the form 

hnZ.^T^hiC/T), (1) 
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where ( represents the set of the c chemical potentials conjugate to the U(l) charges of the 
fluid, V and T represent the volume and the overall temperature of the fluid respectively 
and d — D — 1 is the spacetime dimensions of the boundary. As we have explained above, 
the as yet unknown function h((/T) may be read off from the thermodynamics of large, 
charged, static black holes in AdS. 

The thermodynamic equation of state described above forms an input into the relativis- 
tic Navier-Stokes equations that govern the effective dynamics of the boundary conformal 
fluid. The full equations of fluid dynamics require more data than just the equation of 
state; for example we need to input dissipative parameters like viscosity. However, fluid 
dynamics on S^ 1 admits a distinguished c + n + 1 parameter set of stationary solutions 
(the parameters are their energy E, c commuting charges Ri and n = iank(SO(d)) = [|] 
commuting angular momenta^] on S^ -1 ). These solutions are simply the configurations into 
which any fluid initial state eventually settles down in equilibrium. They have the feature 
that their form and properties are independent of the values of dissipative parameters. 

Although these solutions are nonlinear (i.e. they cannot be thought of as a small fluctu- 
ation about a uniform fluid configuration), it turns out that they are simple enough to be 
determined explicitly. These solutions turn out to be universal (i.e. they are independent 
of the detailed form of the function h((/T)). Their thermodynamics is incredibly simple; 
it is summarised by the partition function 



In Z gc = In Tr exp 



(H — QRi — £l a L a , 



V d T d -'h((/T) 

n: =1 (i-^r u 



where H,L a and Q a represent the energy, angular momenta and the angular velocities of 
the fluid respectively and Vd = Vol(S' a! ~ 1 ) is the volume of the sphere S l . 

We now turn to the gravitational dual interpretation of the fluid dynamical solutions 
we have described above. A theory of a rank c gauge field, interacting with gravity on 
AdS/?, possesses a c + n + 1 parameter set of black hole solutions, labelled by the conserved 
charges described above. We propose that these black holes (when large) are dual to the 
solutions of fluid dynamics described in the previous paragraph. Our proposal yields an 
immediate prediction about the thermodynamics of large rotating black holes: the grand 
canonical partition function of these black holes must take the form of (j2J). 

Notice that while the dependence of the partition function ([2]) on Q/T is arbitrary, its 
dependence on Q a is completely fixed. Thus, while our approach cannot predict thermody- 
namic properties of the static black holes, it does allow us to predict the thermodynamics 
of large rotating black holes in terms of the thermodynamics of their static counterparts!! 
Further, our solution of fluid dynamics yields a detailed prediction for the boundary stress 
tensor and the local charge distribution of the corresponding black hole solution, which 



1 Here, we use the notation [x] to denote the integer part of the real number x. For a list of notation, 
see appendix iDl 

2 The analogue of our procedure in an asymptotically flat space (which we unfortunately do not have) 
would be a method to deduce the thermodynamics and other properties of the charged Kerr black hole, 
given the solution of static charged black holes. 
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may be compared with the boundary stress tensor and currents calculated from black hole 
solutions (after subtracting the appropriate counterterms ji-10]). 

Our proposal is highly reminiscent of the membrane paradigm in black hole physics 



(see, for instance, [Hi, 1121] )■ However, we emphasise that our fluid dynamical description of 
black holes is not a guess; our proposal follows directly from the AdS / CFT correspondence 
in a precisely understood regime (see 13] for a review of AdS/CFT correspondence) |f] 



We have tested the thermodynamical predictions described above on every class of 
black hole solutions in AdS d spaces that we are aware of. These solutions include the most 



general uncharged rotating black holes in arbitrary AdS d space |14H16| , various classes of 



charged rotating black holes in AdS§ x S 5 17H20I]. in AdS-j x S 4 and in AdS^ x S 7 21 



231 ] . In the strict fluid dynamical limit, the thermodynamics of each of these black holes 
exactly reproduces^ ([2]). In all the cases we have checked, the boundary stress tensor and 
the charge densities of these black holes (read off from the black hole solutions using the 
AdS/CFT dictionary) are also in perfect agreement with our fluid dynamical solutions. 
The agreement described in this paragraph occurs only when one would expect it to, as 
we now explain in detail. 

Recall that the equations of fluid mechanics describe the evolution of local energy den- 
sities, charge densities and fluid velocities as functions of spatial position. These equations 
are applicable only under certain conditions. First, the fluctuations about mean values (of 
variables like the local energy density) must be negligible. In the situations under study 
in our paper, the neglect of fluctuations is well justified by the 'large iV' limit of the field 
theory, dual to the classical limit of the gravitational bulk. 

Second, the equations of fluid mechanics assume that the fluid is in local thermodynamic 
equilibrium at each point in space, even though the energy and the charge densities of the 
fluid may vary in space. Fluid mechanics applies only when the length scales of variation 
of thermodynamic variables - and the length scale of curvatures of the manifold on which 
the fluid propagates - are large compared to the equilibration length scale of the fluid, a 
distance we shall refer to as the mean free path l m f p . 

The mean free path for any fluid may be estimated as |24j Z m f p ~ 2 where rj is the shear 
viscosity and p is the energy density of the fluid. For fluids described by a gravitational 
dual, rj = j- where s is the entropy density [241 ] . Consequently, for the fluids under study 



in this paper, / mfp ~ As we will see in § §3.7[ the most stringent bound on / mfp , for 
the solutions presented in this paper, comes from requiring that / m f p be small compared to 
the radius of the S 8 ' 1 , which we set to unity. Consequently, fluid dynamics should be an 
accurate description for our solutions whenever j^- C 1. In every case we have studied, 
it turns out that this condition is met whenever the horizon radius, Rh, of the dual black 
hole is large compared to the AdS radius , i?Ads- Black holes that obey this condition 
include all black holes whose temperature is large compared to unity, but also includes 



3 Alternatively, one could regard the agreement between fluid dynamics and gravity described below as 
a test of the AdS / CFT correspondence (provided we are ready to assume in addition the applicability of 
fluid mechanics to quantum field theories at high density). 

4 See, however, ^6.10l for a puzzle regarding the first subleading corrections for a class of black holes in 
AdS 5 . 
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large radius extremal black holes in AdS§ x S 5 , AdS-j x S 4 and AdS± x S 7 . It, however, 
never includes supersymmetric black holes in the same backgrounds, whose horizon radii 
always turn out to be of the same order as the AdS radius. 

It follows that we should expect the Navier-Stokes equations to reproduce the thermo- 
dynamics of only large black holes. In all the cases we have studied, this is indeed the 
case. It is possible to expand the formulae of black hole thermodynamics (and the stress 
tensor and charge distribution) in a power series in R^ dS /R H . While the leading order 
term in this expansion matches the results of fluid dynamics, we find deviations from the 
predictions of Navier-Stokes equations at subleading orders. 

Starting with the work of Policastro, Son and Starinets [25|], there have been several 
fascinating studies over the last few years, that have computed fluid dynamical dissipative 
and transport coefficients from gravity (see the review [24| and the references therein). 
The work reported in this paper differs from these analyses in several ways. Firstly, the 
solutions of fluid mechanics we study are nonlinear; in general they cannot be thought of as 
small fluctuations about the uniform fluid configuration dual to static black holes. Second, 
all our solutions are stationary: dissipative parameters play no role in our work. 

Indeed our work rather follows the same line of investigation as the one applied to 



plasmaballs and plasmarings in [26|, |27j. These investigations used the boundary fluid 
dynamics to make detailed predictions about the nature and phase structure of the black 
holes and black rings in Scherk-Schwarz compactified AdS spaces. The predictions of 
these papers have not yet been quantitatively verified as the corresponding black hole 
and black ring solutions have not yet been constructed. The perfect agreement between 
fluid dynamics and gravity in the simpler and better studied context of this paper lends 
significant additional support to those predictions of [26|, |27| that were made using fluid 
mechanics. 

While in this paper we have used fluid dynamics to make predictions for black hole 
physics, the reverse view point may also prove useful. Existing black hole solutions in 
AdS spaces provide exact equilibrium solutions to the equations of fluid dynamics to all 
orders in Z m f p . A study of the higher order corrections of these solutions (away from the 
^mfp - * limit) might yield useful information about the nature of the fluid dynamical 
approximation of quantum field theories. 

The plan of our paper is as follows - In §21 we set up the basic fluid mechanical frame- 
work necessary for our work. It is followed by §3] in which we consider in detail a specific 
example of rigidly rotating fluid - a conformal fluid in S 3 xl. A straightforward gen- 
eralisation gives us a succinct way of formulating fluid mechanics in spheres of arbitrary 
dimensions in §U 

We proceed then to compare the fluid mechanical predictions with various types of 
black holes in arbitrary dimensions. First, we consider uncharged rotating black holes in 
arbitrary dimensions in SjB Their thermodynamics, stress tensors and charge distributions 
are computed and are shown to exactly match the fluid mechanical predictions. In §6], we 
turn to the large class of rotating black hole solutions in AdS$ x S 5 . Many different black 
holes with different sets of charges and angular momenta are considered in the large horizon 
radius limit and all of them are shown to fit exactly into our proposal in the strict fluid 
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dynamical limit. However we also find deviations from the predictions of the Navier-Stokes 
equations at first subleading order in / mfp for black holes with all SO (6) Cartan charges 
nonzero (these deviations vanish when the angular velocities, or one of the SO (6) charges 
is set to zero). This finding is at odds with naive expectations from fluid dynamics, which 
predict the first deviations from the Navier-Stokes equations to occur at 0(1^^) and is an 
as yet unresolved puzzle. 

This is followed by §0 dealing with large rotating black holes in AdS± x S 7 and AdSy x S* 4 
backgrounds which are dual to field theories on M2 and M5 branes respectively. The 
thermodynamics of the rotating black hole solutions in these spaces are derived from their 
static counterparts using the duality to fluid mechanics and it is shown how the known 
rotating black hole solutions agree with the fluid mechanical predictions in the large horizon 
radius limit. In each of these cases, the formulae of black hole thermodynamics deviate 
from the predictions of the Navier-Stokes equations only at 0(1^^) in accord with general 
expectations. In the final section, we conclude our work and discuss further directions. 

In appendix |AJ we discuss the constraints imposed by conformal invariance on the 
equations of fluid mechanics. In appendix[Bl we discuss the thermodynamics of free theories 
on spheres. In appendix (Cj we present our computations of the boundary stress tensor for 
two classes of black holes in AdS spaces. In appendix [D], we summarise the notation used 
throughout this paper. 



2 The equations of fluid mechanics 

2.1 The equations 

The fundamental variables of fluid dynamics are the local proper energy density p, local 
charge densities rj and fluid velocities = 7(1, v). Assuming local thermodynamic equi- 
librium, the rest frame entropy density s, the pressure V, the local temperature T and 
the chemical potentials /ij of the fluid can be expressed as functions of p and using the 
equation of state and the first law of thermodynamics, as in § §2. 51 § §2.61 below. In what 
follows, we will often find it convenient to express the above thermodynamic quantities in 
terms of T and /ij rather than as functions of p and ?v 

The equations of fluid dynamics are simply a statement of the conservation of the stress 
tensor T^ v and the charge currents Jf . 

v M r^ = d^T^ + t^ x t Xu + r; A r^ A = , 

2.2 Perfect fluid stress tensor 

The dynamics of a fluid is completely specified once the stress tensor and charge currents 
are given as functions of p, T{ and u** '. As we have explained in the introduction, fluid 
mechanics is an effective description at long distances (i.e, it is valid only when the fluid 
variables vary on distance scales that are large compared to the mean free path / m f p ). As 
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a consequence it is natural to expand the stress tensor and charge current in powers of 
derivatives. In this subsection we briefly review the leading (i.e. zeroth) order terms in this 
expansion. 

It is convenient to define a projection tensor 

projects vectors onto the 3 dimensional submanifold orthogonal to w M . In other words, 
P^ u may be thought of as a projector onto spatial coordinates in the rest frame of the 
fluid. In a similar fashion, —u^u u projects vectors onto the time direction in the fluid rest 
frame. 

To zeroth order in the derivative expansion, Lorentz invariance and the correct static 
limit uniquely determine the stress tensor, charge and the entropy currents in terms of the 
thermodynamic variables. We have 

perfect = fm^ + VP^, 
(^perfect = T&lP, (5) 
(^s) perfect = SU^ , 

where p = p(T, /Xj) is the rest frame energy density, s = s(T, //j) is the rest frame entropy 
density of the fluid and /ij are the chemical potentials. It is not difficult to verify that 
in this zero-derivative (or perfect fluid) approximation, the entropy current is conserved. 
Entropy production (associated with dissipation) occurs only at the first subleading order 
in the derivative expansion, as we will see in the next subsection. 

2.3 Dissipation and diffusion 

Now, we proceed to examine the first subleading order in the derivative expansion. In 
the first subleading order, Lorentz invariance and the physical requirement that entropy 
be non-decreasing determine the form of the stress tensor and the current to be (see, for 



example, §§14.1 of 28]) 



Native = -CM* 4 " - W + <fu» + ttV, 



where 



( — n 

\"% /dissipative H; 
\ S /dissipative ^j- 



= VX. 



(6) 



<f = - K P^(d v T + a v T), 
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are the acceleration, expansion, shear tensor, heat flux and diffusion current respectively. 

These equations define a set of new fluid dynamical parameters in addition to those of 
the previous subsection: ( is the bulk viscosity, rj is the shear viscosity, k is the thermal 
conductivity and Dij are the diffusion coefficients. Fourier's law of heat conduction q = 
— kVT has been relativistically modified to 

q» = -KP^(d v T + a v T), (8) 

with an extra term that is related to the inertia of flowing heat. The diffusive contribution 
to the charge current is the relativistic generalisation of Fick's law. 

At this order in the derivative expansion, the entropy current is no longer conserved; 
however, it may be checked that 

TV,r s = ^ + T{D- l fq^ + (d 2 + 2 W ^. (9) 

As q^, and are all spacelike vectors and tensors, the RHS of ([91) is positive provided 
77, (, k and D are positive parameters, a condition we further assume. This establishes that 
(even locally) entropy can only be produced but never destroyed. In equilibrium, V M J^ 
must vanish. It follows that, g M , gf , 9 and <j^ u each individually vanish in equilibrium. 

From the formulae above, we see that the ratio of T^- tive ^° ^perfect * s °f the or der 
T)l (Rp) where r] is the shear viscosity, p is the rest frame energy density and R is the typical 
length scale of the flow under consideration. Consequently, the Navier-Stokes equations 
may be thought of as the first term in a series expansion of the microscopic equations in 
lmip/R > where Z m f p ~ - p . In this sense, Z m f p plays a role analogous to the string scale in 
the derivative expansion of the effective action in string theory. This length scale may 
plausibly be identified with the thermalisation length scale of the fluid! 

When studying fluids on curved manifolds (as we will proceed to do in this paper), 
one could add generally covariant terms, built out of curvatures, to the stress tensor. For 
instance, we could add a term proportional to R^ v to the expression for T^ v . We will 
ignore all such terms in this paper for a reason we now explain. In all the solutions of fluid 
mechanics that we will study, the length scale over which fluid quantities vary is the same 
as the length scale of curvatures of the manifold. Any expression built out of a curvature 
contains at least two spacetime derivatives of the metric; it follows that any contribution 
to the stress tensor proportional to a curvature is effectively at least two orders subleading 
in the derivative expansion, and so is negligible compared to all the other terms we have 
retained in this paper. 

2.4 Conformal fluids 

We will now specialise our discussion to a conformal fluid - the fluid of the 'stuff' of a 
conformal field theory in d dimensions. Conformal invariance imposes restrictions on both 

J This may be demonstrated within the kinetic theory, where ^ m f p is simply the mean free path of 
colliding molecules, but is expected to apply to more generally to any fluid with short range interactions. 



9 



the thermodynamics of the fluid and the derivative expansion of the stress tensor discussed 
in the previous subsection. 

To start with, conformal invariance requires that the stress tensor be traceless0 This 
requirement relates the pressure of a conformal fluid to its density as V = -j^t (this 
requirement may also be deduced from conformal thermodynamics, as we will see in the 
next subsection) where d is the dimension of the spacetime in which the fluid lives. Further, 
the tracelessness of the stress tensor also forces the bulk viscosity, (, to be zero. 

It is easy to verify that these constraints are sufficient to guarantee the conformal invari- 
ance of the fluid dynamical equations listed above. Consider a conformal transformation 
9ia> = e2 ^9fiu under which fluid velocity, temperature, rest energy density, pressure, entropy 
density and the charge densities transform as 



T = e-*T. 



s = e-^S 



It is easy to verify that these transformations induce the following transformations on 
the stress tensors and currents listed in the previous subsection^ 

Jf = e-^jf, (10) 

_ -d</> jH 

J S ~ e J s- 

These are precisely the transformation properties that ensure the conformal invariance of 
the conservation equations ([3]). See appendix [A] for more details. 



2.5 Conformal thermodynamics 

In this subsection, we review the thermodynamics of the conformal fluids we discuss below. 
The notation set up in this subsection will be used through the rest of this paper. 
Define the thermodynamic potential 

$ = s - rs - HiKi . (ii) 

6 More accurately, conformal invariance relates the nonzero trace of the stress tensor to certain curvature 
forms; for example, in two dimension T£ = -p^R where R is the scalar curvature. However, as we have 
described above, curvatures are effectively zero in the one derivative expansion studied in this paper. All 
formulae through the rest of this paper and in the appendices apply only upon neglecting curvatures. We 
thank R. Gopakumar for a discussion of this point. 

7 Note that under such a scaling, the viscosity, conductivity etc. scale as k = e~^ d ~ 2 ^k , 77 = e~( d ~ 1 ^fj, 
m = e-^fii and D t3 = e'^-^Dij. 
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for which the first law of thermodynamics reads 

d$ = -SdT - VdV - Kidni . (12) 

Let us define Ui = jJLi/T. It follows from conformal invariance and extensivity that 

$ = -VT d h[v) , (13) 

where h(u) is defined by this expression. All remaining thermodynamic expressions are 
easily determined in terms of the function h{vj) 

p=(d-l)V = {d-l)T d h{v), 

n = T d - l h i {v), (14) 

s = T d -\dh-u i h i ) , 

where 

dh 

denotes the derivative of h with respect to its i th argument. 
2.6 A thermodynamic identity 

We will now derive a thermodynamic identity that will be useful in our analysis below. 
Define 

T = S -TS + PV - fuTZi =$ + VVi (15) 
the first law of thermodynamics implies that 

dr = -SdT + VdV - TZidfii . (16) 

Consider scaling the system by a factor (1 + e). Under such a scaling, extensivity implies 
that 

dr = er , dT = dV = d^ = , 

which when substituted into (TIB]) tells us that T = 0. Then we can divide (TI5]) and (TIB]) 

by V to get 

p + V = sT + fan, ^ 
dV = sdT + Tidfii . 

3 Equilibrium configurations of rotating conformal flu- 
ids on S 3 

In this section and in the next, we will determine the equilibrium solutions of fluid dynamics 
equations for conformal fluids on spheres of arbitrary dimension. In this section, we work 
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out the fluid dynamics on S 3 plus a time dimension in detail^ In the next section, we 
generalise the results of this section to spheres of arbitrary dimension. 

3.1 Coordinates and conserved charges 

Consider a unit S 3 embedded in IR 4 as 

x 1 = sin 9 cos <j>i 

x 2 = sin 9 sin <b\ , , 

a \ (18) 

X = COS 9 COS <p2 

x 4 = cos 9 sin <p 2 

with 9 G [0, |], 4> a G [0, 27r). The metric of the spacetime S* 3 x IR is 

ds 2 = -dt 2 + d9 2 + sin 2 9 d(f)j + cos 2 9 d<j) 2 2 . (19) 

This gives the following non-zero Christoffel symbols: 

rji* = ~ r 5 2 2 = - cos ^ sin 61 , r* = T% e = cot 0, rfo = r£, = - tan . (20) 

For the stationary, axially symmetric configurations under consideration, dtT^ = 
d^T^ = 0. Using p|), © becomes 

= V M T^ = d e T et + (cot - tan 9)T et , (21) 

= V M T^ = <9 e T ee + (cot 9 - tan 0)T* fl + cos 9 sin (T^ 1 - T^ 2 ) , (22) 

= V^T^ 1 = deT 9 ^ 1 + (cot - tan 0)T^ + 2 cot T^ 1 , (23) 

= V M T^ 2 = <9 e T^ 2 + (cote - tan#)T e<fc - 2tan#T e<fc . (24) 

The Killing vectors of interest are d t (Energy) and d^ a (SO(4) Cartan angular mo- 
menta). Using the formula for the related conserved charge, J d 3 x^/—g T°^g fiv k 1 ' ', we get: 

E = J dOdfadfc cos 9 sin 9T a , 

L x = JdBdfadfa cos 9 sin 3 9 T^ 1 , (25) 

L 2 = J dBdfadfc cos 3 9 sin 9 T^ 2 . 

Assuming = = (as will be true for stationary solutions we study in this paper), the 
entropy and the R-charges corresponding to the currents in ([5]) are given as 



S = I d9d(fiid(j)2 cos 9 sin 9 7 s , 

f (26) 
Ri = d6*d0id02 cos 9 sin 3 9 77^ . 



8 In this case, the dimensions of the spacetime in which the fluid lives is d = 3 + 1 = 4. The number 
of mutually commuting angular momenta is n = 2. The black hole dual lives in AdS space of dimensions 
D = d+ 1 = 5. 
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3.2 Equilibrium solutions 



As we have explained in the § §2. 31 each of the three quantities o^ u , g M , gf must vanish on 
any stationary solution of fluid dynamics. The requirement that = has a unique 
solution - the fluid motion should be just a rigid rotation. By an 50(4) rotation we 
can choose the two orthogonal two planes of this rotation as the (1-2) and (3-4) planes 
(see ffl8]) ). This implies that m m = 7(1, 0, co>i, co> 2 ) (where we have listed the (t, 9, (fix, fa) 
components of the velocity) with 7 = (1 — v 2 ) and v 2 = <J\ sin 2 9 + io\ cos 2 9, for some 
constants u>i and uj 2 . 

Our equilibrium fluid flow enjoys a symmetry under translations of t, 0i and fa; con- 
sequently all thermodynamic quantities are functions only of the coordinate 9. 
Evaluating the tensors in ([7j), we find 

a" = (0, -ft In 7, 0,0), 
■& = 0, 
a*" = 0, 



r 



7 

T 



,0,0 
,0,0 



(27) 



The requirement that and vanish forces us to set 

T = T7 , fii = Tvi, 



(2f 



for constant r and z/j. These conditions completely determine all the thermodynamic 
quantities as a function of the coordinates on the sphere. We will now demonstrate that 
this configuration solves the Navier-Stokes equations. 

First note that for an arbitrary rigid rotation, the dissipative part of the stress tensor 
evaluates to 



rpjJLV 

dissipative 



-«7 



/o 
1 



1 









ui 
\0 u 2 












d 


"T" 


d9 


.7. 



(29) 



an expression which simply vanishes once we impose (128]) . Consequently, all nonzero con- 
tributions to the stress tensor come from the 'perfect fluid piece' and are given by 



T^ v = -y 2 

perfect " 



f(p + v 2 V) oj^p + V) oo 2 (p + V) \ 

7~ 2 P 

wi (p + V) Ljfp + (esc 2 9-u)\ cot 2 9)V lo x lo 2 (p + V) 

\uj 2 {p + V) uiu 2 {p + V) u\p + (sec 2 9 - cuf tan 2 9)V J 

(30) 
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The only non-trivial equation of motion, (122]) . can be written as 



dV p + Vd-f 



d6 



7 



d9 



0. 



(31) 



Now using the thermodynamic identity (fl71) we may recast (I3T]) as 



7S d^ 



T 

7. 



+ in d9 



/i, 

7 



0. 



(32) 



an equation which is automatically true from ( 128]) . Consequently, rigidly rotating config- 
urations that obey (128]) automatically obey the Navier-Stokes equations. 

In a similar fashion, it is easy to verify that all nonzero contributions to the charge 
currents come from the perfect fluid piece of that current, and that the conservation of 
these currents holds for our solutions. 

In summary the 3 + c parameter set of stationary solutions to fluid mechanics listed 
in this subsection (the parameters are r, u a and z/j where i — 1 . . . c ) constitute the most 
general stationary solutions of fluid mechanics. 



3.3 Stress tensor and currents 

Using the equations of state ( lT4"j) . we find that 

p = 3V = 3r 4 7 4 /i(^), 



s = T 3 7 3 [4/i(z/) - Uihiiu)}, 
n = r 3 7 3 /ii(z/). 



(33) 



The stress tensor is 

/3 + v 2 

= r 4 A 7 6 



4ui 

1 - v 2 
Alox 3uo 2 + esc 2 9 - ujI cot 2 9 



4o; 2 




\ 



\ 4o> 2 AujiUJ 2 
Charge and entropy currents are given by 



• (34) 



3o;| + sec 2 9 - tuf tan 2 9 J 



J? = r 3 7 4 a(l,(W, W2 ) 



J^rVE^O,^,^), 



where we have defined 



A = h{v) , 

B = Ah(u) - Vihi{v) , 
Ci = hi(v) = — . 



(35) 



(36) 
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3.4 Charges 



The energy, angular momentum, entropy and R-charges may now easily be evaluated by 
integration: we find 

2ul 



E 



Li 



V4T A 
V A r A A 



+ 3 



V 4 t 4 A 



2wi 



S 



2u; 2 



(37) 



w ?)(i 



where V4 = Vo^S 13 ) = 27r 2 is the volume of S" 3 . These formulae constitute a complete 
specification of the thermodynamics of stationary rotating conformal fluids on S 3 . 



3.5 Potentials 

In the previous subsection we have evaluated all the thermodynamic charges of our rotat- 
ing fluid solutions. It is also useful to evaluate the chemical potentials corresponding to 
these solutions. To be specific we define these chemical potentials via the grand canonical 
partition function defined in the introduction 

Z gc = Trexp (±(-H + Q a L a + ( t R t )^ = exp (- E - TS ~^ aLa ~ ^ , (38) 

where the last expression applies in the thermodynamic limit. In other words 

r= ©W "^Ow Ci= O s , il , K / (39) 

It is easy to verify tha1@ 

T = t, Q a = u a , C« — TV % ■ (40) 

Note that T, Q a and Q are distinct from T, u a and /ij. While the former quantities are 
thermodynamic properties of the whole fluid configuration, the latter quantities are local 
thermodynamic properties of the fluid that vary on the S 3 . In a similar fashion, the energy 
E of the solution is, of course, a distinct concept from the local rest frame energy density p 
which is a function on the sphere. In particular, E receives contributions from the kinetic 
energy of the fluid as well as its internal energy, E. 

9 We can express dE — rdS — uj a dL a — TVidRi in terms of dr, doj a , dvi and check that it vanishes, or we 
can check the Legendre transformed statement d{E — tS — Lu a L a — TViR{) = —Sdr — L a daj a — Rid(ri>i). 
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3.6 Grand canonical partition function 




(41) 



c (1-0?)(1-01) 



where V4 = V(S 3 ) = 2tt 2 is the volume of S 3 . 

In other words, the grand canonical partition function of the rotating fluid is obtained 
merely by multiplying the same object for the non-rotating fluid by a universal angular 
velocity dependent factor. 

3.7 Validity of fluid mechanics 

A systematic way to estimate the domain of validity of the Navier-Stokes equations would 
be to list all possible higher order corrections to these equations, and to check under what 
circumstances the contributions of these correction terms to the stress tensor and currents 
are small compared to the terms we have retained. Rather than carrying out such a detailed 
(and worthwhile) exercise, we present in this section a heuristic physical estimate of the 
domain of validity of fluid dynamics. 

Consider a fluid composed of a collection of interacting 'quasiparticles', that move at an 
average speed v p and whose collisions are separated (on the average) by the distance Z m f p 
in the fluid rest frame. Consider a particular quasiparticle that undergoes two successive 
collisions: the first at the coordinate location x± and subsequently at x 2 . In order for the 
fluid approximation to hold, it must be that 

1. The fractional changes in thermodynamic quantities between the two collision points 
(e.g. [T(xi) — T(x2)]/T(xi)) are small. This condition is necessary in order for us to 
assume local thermal equilibrium. 

2. The distance between the two successive collisions is small compared to the curva- 
ture/compactification scales of the manifold on which the fluid propagates. This 
approximation is necessary, for example, in order to justify the neglect of curvature 
corrections to the Navier-Stokes equations. 

Let us now see when these two conditions are obeyed on our solutions. Recall that 
the local temperatures in our solutions take the form T = T7 where T is the overall 
temperature of the solution. If we treat the free path Z m f p as a function of temperature and 
chemical potentials, conformal invariance implies that 



7 

Hence, the first condition listed above is satisfied when the fractional variation in (say) 
the temperature is small over the rest frame mean free path Z m f p (T, z/j), i.e. provided 





(42) 
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which must hold for all points of the sphered! The strictest condition one obtains from 
this is 




V 1 - u l - V 1 -^ 



It turns out that the second condition listed above is always more stringent, espe- 
cially when applied to fluid quasiparticles whose rest frame motion between two colli- 
sions is in the same direction as the local fluid velocity It follows from the formulae of 
Lorentz transformations that the distance on the sphere between two such collisions is 
lmfp(T, ^1)7(1 + v/v p ) = l m { p (T, Vi)(l +v/v p ), where v p is the quasiparticle's velocity in the 
rest frame of the fluid and v the fluid velocity. As the factor (1 + v/v p ) is bounded between 
1 and 2, we conclude that the successive collisions happen at distances small compared to 
the radius of the sphere provided 



Hence, we conclude that the condition (j44H (which is always more stringent than f|43]) ) 
is the condition for the applicability of the equations of fluid mechanics. 

Of course the model (of interacting quasiparticles) that we have used to obtain (jUJ) 
need not apply to the situations of our interest. However the arguments that led to (134"]) 
were essentially kinematical which leads us to believe that the result will be universal. 
Nonetheless, it would be useful to verify this result by performing the detailed analysis 
alluded to at the beginning of this subsection. 

4 Rotating fluids on spheres of arbitrary dimension 

We now generalise the discussion of the previous section to the study of conformal fluids 
on spheres of arbitrary dimension. 
Let us embed S 2n in R 2n+1 as 



Where 9 n G [0, 7r], all other 6 a G [0, |] and 4> a G [0, 2tt). Any products with the upper limit 
smaller than the lower limit should be set to one. Although we appear to have specialised 

10 Recall that all variations in the temperature are perpendicular to fluid velocities, so that the typical 
scale of variation in both the rest frame and the lab frame coincide. 



lmfp(T, Vi) < 1 . 



(44) 




(45) 
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to even dimensional spheres above, we can obtain all odd dimensional sphere, S 
by setting 8 n = ir/2 in all the formulae of this section. 
The metric on S 2n x time is given by 



2n-l 



n I a— I 



n / a— I 



ds 2 = _ dt 2 + cos 2 e b Uel + II cos2 e » ) sin2 e ^ 

o=l \b=l / 0=1 \6=1 

We choose a rigidly rotating velocity 

u l = 7 u 6a = 7/ = 7Co> a 

7 = (1- W 2 )-V2 „ 



n / a— 1 



a=l \6=1 

As in § §3.2[ the equations of motion are solved, without dissipation, by setting 



T 

7 

which gives the densities 



r = constant, 



T 



Vi = constant, 



p=(d-l)P = (d- l)rV^(^), 
r i = r*- 1 7 d - 1 /i i (i/), 

This gives a stress tensor 

T u = r d A{d-f d+2 - l d ) T t<t}a = T*»* = r d AdY +2 uJ a 

/a-l 

T e a e a = r d Al d j -Q gec 2 Qh 



vb=l 



and currents 
where 



'a-l 



di» 2 e£ + Y J] sec 2 6/ 6 esc 2 d 



T <t>a<t> b = T d Ad~i d+2 UJ a UJ b 



J l s = r^By" J s a = Jf° = r^'BYuJa 



d-1 D-,d, 



j* = T^crf J- a = o Jf a = T d ~ x Ca a uj a 

A = h{v) , 

Z? = dh{v) — Vihi{y) , 
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Integrating these giver 11 
E 
S 

L a 

Hi 



V d T d h{v) 

ryi-^ 2 ) 



+ d-l 



V d r d - l [dh(u) - uMu)] 



n 6 (i-w 5 



b ) 

2uj a 



V d T d h(v) 

n 6 (i - ^ ' 



where 



V d = Vol(5 



2-7T d / 2 



Differentiating these gives 

T = r 

and the grand partition function 

InZ, 



r(d/2) ' 
IL(i-ng) ' 



(53) 



(54) 



(55) 



As in the previous subsection, the fluid dynamical approximation is expected to be valid 
provided Z mfp (T, Ui) < 1. 

In appendix [Bj we have computed the thermodynamics of a free charged scalar field on 
a sphere, and compared with the general results of this section. 



5 Comparison with uncharged black holes in arbi- 
trary dimensions 

In the rest of this paper, we will compare the predictions from fluid dynamics derived 
above with the thermodynamics, stress tensors and charge distributions of various classes 
of large rotating black hole solutions in AdS spaces. We start with uncharged rotating 
black holes on D dimensional AdS spaces (where D is arbitrary), which are dual to rotating 
configurations of uncharged fluids on spheres of dimension (D — 2). 

11 In deriving these formulae we have 'conjectured' that J S d-i 7 d = — jy- ^ * s eas y to derive this 

formula for odd spheres. We have also analytically checked this formula for S 2 and S 4 . We are ashamed, 
however, to admit that we have not yet found an analytic derivation of this integral for general even 
spheres. 
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5.1 Thermodynamics and stress tensor from fluid mechanics 



In case of uncharged fluids the function h(u) in the above section is a constant h(u) = h. 

qual to 5 

V d T d h 



Therefore hi{v) = are all equal to zero. It follows from equations (j55|) and that 



E 

S 



n 6 a-^) 

V d T d - l hd 

IL(i-ng)' 



^ 1 - fi2 + 



(56) 



a n 6 (i - « 

R; = 0. 



2a, 



i - a 



The partition function is given by 



InZ 



V d T d ~ l h 



The stress tensor becomes 

T tt = hT d (dY +2 - l d ) T t(/,a = T^ at = hT d dY +2 tt a 



(57) 



T e a e a = hT d ld -Q gec 2 9h 



v6=l 



(5f 



'O-I 



v 6=l 



T 0a0 6 = hT d d~f d+2 n a n b . 



The mean free path in fluid dynamics can be estimated by taking the ratio of shear 
viscosity to energy density. As mentioned in the introduction, for fluids with gravity duals 
we can equivalently estimate / m f p by taking the ratio of entropy to 4tt times the energy 
(because of the universal relation s = Attt]). 



lm{p{T, v) 



S 



AttE 



AiiT{d- 1)' 



(59) 



Consequently the expansion in Z m f p translates simply to an expansion in inverse powers of 
the temperature of our solutions. 

5.2 Thermodynamics from black holes 

The most general solution for uncharged rotating black holes in AdSz) was obtained in 
15 . 16]. These solutions are labelled by the n + 1 parameter^ ctj and r + (these are related 



12 Recall that n denotes the number of commuting angular momenta and is given by the expression 
n = rank [SO(D - 1)] on AdS D . 



20 



to the n angular velocities and the horizon radius (or equivalently the mass parameter) of 
the black holes). The surface gravity k and the horizon area A of these black holes are 
given b} 



r+(l + rl 



r+(l + rl 



n 1 



1 



i=l 



r+ + a\ r + 

1-ri 



^.2 i_ 



i=l 



ri + a- 



2r, 



when D = 2n + 1 , 



when D = 2n + 2 



.4 



i=l 



ri + a' 



(60) 



-a? 



n 2 i 2 



mi 



when Z) = 2n + 1 , 
when D = 2n + 2 . 



i=l 



We will be interested in these formulae in the limit of large r + . In this limit the parameter 
m (which appears in the formulae of 15, [3]) and the temperature T = k/2tc are given as 
functions of r + by 



T 



47T 



(l + 0(l/^)) 



(61) 



2m = r°~ 1 (l + 0(l/r*)) . 
From these equations, it follows that the parameter m is related to the temperature T as 

47T 



2m 



D-l 



(1 + 0(1/T 2 )). 



D - 1 

To leading order in r + , the thermodynamic formulae take the form 

V D T" 1 r 1 l/ - , - L 



(62) 



Vb-iT 



D-l 



5 



i6vrG D n;=i(i-«p 

y D -iT g ~ 2 (D - 1) 

i67rG D n;=i(i-«^ 




~l D-l 



St 

.i=l 

2 a,; 



2a? 



-a; 



+ D-2 



a? 



(63) 



= 0, 

where Vd_i is the volume of S D ~ 2 and Ge> is Newton's constant in D dimensions. The 
corrections to each of these expressions are suppressed by factors of 0(l/r+) = 0(1/T 2 ) 



13 In the expression of k for even dimension, the sign inside the second term in equation (4.7) of [16[ is 
different form the sign given in equation (4.18) of [15j; we believe the latter sign is the correct one. 
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relative to the leading order results presented above (i.e. there are no next to leading order 
corrections) . 

These thermodynamic formulae listed in ( l63l) are in perfect agreement with the fluid 
mechanics expressions in (1561) upon making the following identifications: the spacetime 
dimensions of the boundary theory d = D — 1, the black hole angular velocities ai are 
identified with Q a and the constant h is identified as 



h 



16nG 



D 



47T -> D - 1 



D - 1 



(64) 



In the next subsection, we will see that this agreement goes beyond the global thermody- 
namic quantities. Local conserved currents are also in perfect agreement with the black 
hole physics. 



5.3 Stress tensor from rotating black holes in AdS 



D 



The uncharged rotating black holes both in odd dimensions (D = In + 1) and even dimen- 



sions (D = 2n + 2) are presented in detail in [15J, equation (E-3) and (16J, equation (4.2 



After performing some coordinate transformations (see appendix 10.11) that take the metric 
of that paper to the standard form of AdS b at the boundary, we have computed the stress 
tensor of this solution. 

Our calculation, presented in appendix [CJ] uses the standard AdS/CFT dictionary. In 
more detail, we foliate the solution in boundary spheres, compute the extrinsic curvature 
0(J of these foliations near the boundary, subtract off the appropriate counter terms con- 
tributions (3 -10], and finally multiply the answer by the r D ~ l to obtain the stress tensor 



on a unit sphere. 

We find that the stress tensor so calculated takes the form (see appendix IC.ll) 

n " = n " = ^ D - l ^ (65, 




Here 7 2 = 1 - J2a=i ^l^l where /i a = (IILi cos ®b) sin 6 a . 

Note that the functional form of these expressions (i.e. dependence of various com- 
ponents of the stress tensor on the coordinates of the sphere) agrees exactly with the 
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predictions of fluid dynamics even at finite values of r + . In the large r + limit (using (|62|) 
and flM|) ) , we further have 



D — 1 = d, 



2m T d h. 



16nG D 



With these identifications, (1651) coming from gravity agrees precisely with (1581) from fluid 
mechanics. 

We proceed now to estimate the limits of validity of our analysis above. From the black 
hole side, since we have expanded the formulae of black hole thermodynamics in l/r + 
to match them with fluid mechanics, this analysis is valid if r + is large. From the fluid 
mechanics side, we expect corrections of the order of Z m f p . To estimate / m f p in this case, we 
substitute (ISTI) into (I59I) to get 



Hence, we see that the condition from fluid mechanics is exactly the same as taking large 
horizon radius limit: the expansion of black hole thermodynamics in a power series in — 
appears to be exactly dual to the fluid mechanical expansion as a power series in Z mfp . 



6 Comparison with black holes in AdS$ x S 5 

Large N, Af — 4 Yang-Mills, at strong 't Hooft coupling on S 3 x R, is dual to classical 
gravity on AdS^ x S 5 . Hence, we can specialise the general fluid dynamical analysis pre- 
sented above to the study of equilibrium configurations of the rotating Af = 4 plasma on 
S 3 and then compare the results with the physics of classical black holes in AdS$ x S 5 . 

Large black holes in AdS§ x S 5 are expected to appear in a six parameter family, labelled 
by three SO (6) Cartan charges (c = 3), two 5*0(4) rotations (n = 2) and the mass. While 
the most general black hole in AdS$ x S 5 has not yet been constructed, several sub-families 
of these black holes have been determined. 

In this section, we will compare the thermodynamic predictions of fluid mechanics with 
all black hole solutions that we are aware of and demonstrate that the two descriptions 
agree in the large horizon radius limit. For one class of black holes we will also compare 
black hole stress tensor and charge distributions with that of the fluid mechanics and once 
again find perfect agreement (in the appropriate limit). 

We begin this section with a review of the predictions of fluid mechanics for strongly 
coupled Af = 4 Yang-Mills on S 3 . Note that this is a special case of the conformal fluid 
dynamics of previous sections with d = D — 1 = 4. 
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6.1 The strongly coupled J\f = 4 Yang-Mills Plasma 

The gravity solution for SO (6) charged black branes (or, equivalently, large SO (6) charged 
but non- rotating black holes in AdS<$ x S 5 ) has been used to extract the equation of state 
of TV = 4 Yang- Mills (see 0, §2] for the thermodynamic expressions in the infinite radius 
limit). 

Rather than listing all the thermodynamic variables, we use the earlier parametrisation 
of (fT4"j) to state our results. The thermodynamics of the M = 4 Yang-Mills is described by 
the following equations^ 

h(v) = —— = 2n N =^ — — , 

# 27rn(l + ^) ^ 



T (2 + Ej «i - n, «*) 1 + * ' ( 66 ) 



W = ^3 



2vriV 2 n,(l + ^) 2 



2-E 7 «i-lI,«. "" 



where the auxiliary parameters Ks have a direct physical interpretation in terms of entropy 
and charge densities (see §2 of |29|) - 



4tt 2 ^ 2 

constrained by > and by the condition^ 

2 + Ej «j - rij 



(67) 



n,(i+«. 



> 0. 



It follows from (1671) that is finite for configurations with finite charge and non-zero 
entropy. The configurations with — > oo (for any 2) are thermodynamically singular, 
since in this limit, the i th charge density is much larger than the entropy density. Hence, 
in the following, we shall demand that Ki be finite. 

The general analysis presented before now allows us to construct the most general 
stationary solution of the M = 4 fluid rotating on a 3-sphere. The thermodynamic formulae 
and currents of these solutions follow from (j55|) . (1541) and (1571) upon setting 

H(i 1 ! 



A = h{u) = 2n z N' 



(2+E ? s-n ? s 



B = 4 % ) - = ^ T^W^TO ' ^ 



■(2+E i «i-n i « 



3 



14 Note that our convention for the gauge field differs from [23, §2] by a factor of \/2. 

15 Which is obtained by requiring that the temperature T > in the expression for T in [^J, §2]. 
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which leads to 



and 



2ttT n ■(! + «; 



InZ, 



2vr 2 iV 2 \/4T 3 n j (l + ^) 3 



i • 



where we have used the notation V4 = Vo^S* 3 ) = 2n 2 as before. 

As before, the mean free path in fluid mechanics can be estimated as 



'■mfp 



s 



AtiE 



B 



(2 + J2j *i - Yij 



n=0 {d-l)^TA 67rrn,(i + ^; 
- > - 

+ n] 



QttT 



j 



(69) 



(70) 



(71) 



6.2 The extremal limit 

The strongly coupled M = 4 Yang-Mills plasma has an interesting feature; it has interesting 
and nontrivial thermodynamics even at zero temperature. In this subsection, we investigate 
this feature and point out that it implies the existence of interesting zero temperature 
solutions of fluid dynamics which will turn out to be dual to large, extremal black holes. 



6.2.1 Thermodynamics 

In the above section, we presented the thermodynamics of strongly coupled M = 4 Yang- 
Mills plasma in terms of the parameters K{. These parameters are constrained by the 
conditions Ki > and y+^- > 1 with Ki finite. In order to visualise the allowed range 
over which the variables fi^'s can vary, it is convenient to define a new set of variables 

Xi = , Xi = X, Y, Z, 

* ji (72) 

X= X + Y + Z - 1 ' 

The constraints Ki > and Y] 4 tt— > 1 with K; finite translate into the constraints 
< Xi < 1 and X + Y + Z > 1 . Geometrically, this is just the statement that Xj's can 
lie anywhere inside the cube shown in figiH away from the planes Xj = and on or above 
the plane X + Y + Z = 1. 
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The energy density, the entropy density and the charge densities of the Yang-Mills 
plasma may be rewritten as a function of X, Y, Z and x as 



p = £m 2 N 2 XYZ X *, s = An z N' z XYZx 



2 AT 2- 



n = 27iN 2 XYZ X i 



1-Xi 



(73) 



Xi 



The condition for the validity of fluid mechanics becomes 

1 



mfp 



6ttx 



< 1 or x > 1 • 



(74) 




Figure 1: The space of allowed /%'s. The axes correspond to X 



l+Kl 



Y 



1 + K2 



and Z 



1+K3 



. The Xj's can lie anywhere in the cube outside the "extremal" plane X + Y + Z = 1. 



Consider now the case in which X is large, but finite and X, Y, Z take values close to 
the interior of the triangle ABC in fig {H From (1721) and f!73|) . it is evident that this is 
equivalent to taking an extremal limit T —>■ with appropriate chemical potentials. All 
thermodynamic quantities listed above are smooth in this limit and the fluid mechanics 
continues to be valid. 

The M = 4 Yang-Mills plasma with three nonzero R-charges always has a nonsingular 
extremal limit. In the case that one of the charges say r% is zero, then we are constrained 
to move on the X% = 1 plane in the space of Xi's. Hence, we can never approach the 
'extremal triangle' X + Y + Z = l@ Thus, we have no nonsingular extremal limit if any 
one of the three R-charges is zero. By a similar argument, no nonsingular extremal limit 
exists if two of the R-charges were zero. 

We note that Gubser and Mitra have previously observed that charged black branes 

Se- 



near extremality are sometimes thermodynamically unstable 



Although we have not 



16 Remember that we have already excluded, on physical grounds, the point X\ 
lies in the intersection of X$ = 1 plane and the extremal plane X + Y + Z = 1. 



x 2 = 0, x 3 



1 which 
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performed a careful analysis of the thermodynamic stability of the charged fluids we study 
in this paper (see however [iij]), we suspect that these fluids all have Gubser-Mitra type 
thermodynamic instabilities near extremality. If this is the case, the near extremal fluid 
solutions we study in this section and the next - and the black holes that these are dual to - 
are presumably unstable to small fluctuations. Whether stable or not, these configurations 
are valid solutions of fluid dynamics. We postpone a serious discussion of stability to future 
work0 



6.2.2 Fluid mechanics 



The thermodynamic expressions for the charges of a rotating Yang-Mills plasma take the 
form 
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2tv 2 N 2 XYZVa 
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(75) 



2wN 2 XYZV d 



(l-atf)(l 
and the mean free path 
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X, 
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mfp 



6ttT 



<C 1 



(76) 



We see that all thermodynamical charges of our rotating fluid configurations are nonsingu- 
lar, and that fluid mechanics is a valid approximation for these solutions, in the extremal 
limit described in the previous subsection, provided only that x ^ 1@ 

The solution so obtained describes a rotating fluid whose local temperature vanishes 
everywhere, but whose rest frame charge density is a function of location on the S 3 (it 
scales like 7 3 ). As we will see below these extremal configurations of rotating fluid on S 3 
are exactly dual to large, rotating, extremal black holes in AdSs. 



6.3 Predictions from fluid mechanics in special cases 

As mentioned in the beginning of this section, the most general black hole in AdS^ x S 5 
has not yet been constructed, but several subfamilies of these black holes are known. To 

17 We thank Sangmin Lee for discussion of these issues. 

18 In greater generality, in order for fluid mechanics to be a valid approximation for our solutions it is 
necessary that either T ^> 1 (which is by itself sufficient) or that X + Y + Z — 1 — > (under which condition 
the ratio x °f the previous section must be large and (conservatively) none of X, Y or Z be very small). 
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facilitate the comparison between fluid mechanics on S 3 on one hand and these subfamilies 
of black holes on the other, in this subsection, we specialise the general predictions of the 
previous subsection to various specific cases. 



6.3.1 All SO (6) charges equal: arbitrary angular velocities 

Consider first the case of a fluid with equal SO (6) charges (with the rotational parameters 
arbitrary). That is we set K\ = k,2 = ^3 = k in the general formulae above. Noting that 
(2 + 3k — k 3 ) — (k + 1) 2 (2 — k) we find that the stress tensor and currents are given by 
dMD and d35]) with 

„ 27T 2 iV 2 (l + k) 

A = 



(2 - Kf ' 

4tt 2 N 2 

a 2 * N2 



(2 - Kf ■ 

The thermodynamics can be summarised by 

lifTyfR 2^N 2 V A T\l + k) 

0= (2^)' l^ g c(T,fi,C)= (1 _ fi?)(1 _^ )(2 _ K)4 . (78) 

The formula for mean free path (ITTj) reduces to 

2 — K 



^mfp ~ 



1 + K 



(79) 



Let us specialise the extremal thermodynamics of M = 4 Yang-Mills fluid presented 
before to this case. In terms of the variables introduced in § §6. 21 we have X = Y = Z 
which is a straight line in the Xi space. The extremal limit is obtained when this line cuts 
the extremal plane X + Y + Z = 1 , i.e, at the point X — Y — Z — 1/3. This corresponds 
to the extremal limit k —>■ 2. 

More explicitly, in the extremal limit 

T^O, (2-«) = -, (80) 

X 

with x large but finite. The thermodynamic quantities obtained by differentiating the 
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grand partition function (1931) . 

2N 2 n 2 V 4 T 3 



S 
L x 

U 
R 
E 



{2-nf (l-a 2 )(l-6 2 ) 
4,N 2 Ti 2 V A T A {l + k) a 

(2 - k) a (1 -a 2 ) 2 (l -b 2 ) 
AN 2 7T 2 V 4 T 4 (1 +k) b 



(2 - «)* 
27rA^ 2 y 4 r 3 ^ 



l-a 2 )(l-6 2 ) 2 

1 



(2-/t) 3 (l-a 2 )(l-6 2 ) 
27r 2 iV 2 \/4T 4 (l + k) [4 - (1 + a 2 )(l + 6 2 )] 



(2 - k)* 



(1 -a 2 ) 2 (l -6 2 ) 2 



are all smooth; and they describe a fluid configuration whose energy, angular momentum, 
charge and entropy scale as N 2 x A , N 2 x 4 , N 2 x 3 and N 2 x 3 respectively. 

6.3.2 Independent SO(6) charges: equal rotations 

Consider the special case u% = uj 2 = ft (the three SO (6)) chemical potentials are left 
arbitrary). The stress tensor and currents are given by (jMj) and fl35|) with 



A = 2tt 2 N 2 



B = 4n 2 N 2 



n*(i+«i 



(2+E 7 -^-n,«i) 4 



n<(i+*j 



d = 2tcN 2 ^ 
1 

7 



n,(i+«i) 2 



(82) 



The thermodynamics can be summarised by 
2 7 rTn i (l + /S-) 



i + 



InZ. 



2vr 2 iV 2 ^4T 3 n :; (l + ^) 3 



The expression for mean free path (17111 reduces to 



'mfp 



34) 



The extremal limit Ej(l + K j) 1 ~~ > 1 with all Kj kept finite, is nonsingular, and yields 
solutions that are well described by fluid dynamics when / m f p is small. 



29 



6.3.3 Two equal nonzero SO(6) charges: arbitrary angular velocities 

Consider now the case when k\ — K2 — K, K3 — 0. We find that the stress tensor and 
currents are given by ( 134"1) and ( 1351) with 

_ tv 2 N 2 (1 + k) 2 
A ~ 8 ' 



C\ — C*2 — — 

c 3 = o. 

The thermodynamics can be summarised by 



2 ' (85) 
vriV 2 ^! + k) 
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N 2 V 4 T 3 (1 + k) 



Ci = 6 = *zV«, in^(r,n,c)= 8(1 _ si;)( ;_^ ) . (86) 

The expression for mean free path, from (J7T1) . is 

/mfp ~ 37rT(l + «) • (8?) 

It follows from flHTl) that fluid mechanics is a good approximation when T is large. Though 
this equation would appear to suggest that the fluid dynamical approximation is also valid 
(for instance) at fixed T and large k, we have emphasised before, the limit of large k is 
thermodynamically suspect. Conservatively, thus, fluid mechanics applies only at large 
temperatures. 

6.3.4 A single nonzero charge: arbitrary angular velocities 

We now set K\ = k, K2 = ^3 = leaving angular velocities arbitrary. The stress tensor and 
currents are given by (|34"|) and fl35l) with 

2vr 2 iV 2 (i + «;) 3 

(2 + kY ' 
47r 2 iV 2 (l + K ) 2 

(2 + nf ' (88) 
2vriV 2 ^K(l + k) 2 



B 



(2 + Kf 

C 2 = C 3 = . 

The thermodynamics can be summarised by 

r 2nT ^ ]n ? (to n ^ 2 n 2 v,t\i + nf 
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The mean free path, from (ITT!) is given by 



mfp 



QixT 



2 + K 
1 + K 



(90) 



As in the previous subsection, this particular case does not admit thermodynamically 
nonsingular zero temperature (or extremal) configurations. 



6.4 Black holes with all R-charges equal 

Having derived the fluid mechanics predictions for various different black holes, we now 
proceed to examine the black hole solutions. First, we will focus on the case of black holes 
with arbitrary angular momenta in AdSs but equal SO (6) charges. The relevant solution 



has been presented in [18 



6.4.1 Thermodynamics 



The black holes presented in [18] are labelled by two angular velocities a, 6, and three 
more parameters q, m and r + . These five parameters are not all independent; they are 
constrained by one equation relating horizon radius to the parameter m (A r = in that 
paper). We thus have a four parameter set of black holes 1^1 

The relatively complicated black hole thermodynamic formulae of [l8[ simplify if the 
parameter r + (which may be interpreted as the horizon radius) is taken to be large. In 
particular, consider the limit 



^> 1 and y = q/r+ fixed. 



In this limit, to leading order, we have 

T 
2m 



r 4 



(2-y 2 



2tt 

rUi + y 2 



(91) 



(92) 



From the positivity of T and r + it follows immediately that < y 2 < 2. 



Multiplying all thermodynamic integrals in pL8j by = and noting that our 

charge R is equal to their Q/\/3, the black hole thermodynamic formulae reduce to (to 
leading order in r + ) 



a 
b. 



2nyT 



In Za 



(2 

2rt 2 N 2 (l + y 2 ) 



(93) 



V,T 3 



{2-y 2 y L(i-n?)(i-n2)J ' 
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We work in conventions in which the AdS radius and hence the parameter g of [181 ] is set to unity. 
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Once we identify the black hole parameter y 2 with the fluid parameter k, these formula 
take precisely the form of fluid mechanics formulae ([75]) with the equation of state coming 
from (|77J)@ 

We can now compute the fluid mechanical mean free path Z m f p as a function of bulk 
black hole parameters. From equations (1921) and (1791) . we find (assuming that r + is large) 

mfp ~3r + (l + *;)- 

As 1 + k = 1 + y 2 is bounded between 1 and 2, it appears from this equation that the 
expansion in powers of l/r + is simply identical to the fluid dynamical expansion in powers 
of / m f p . This explains why black hole thermodynamics agrees with the predictions of the 
Navier-Stokes equations when (and only when) r + is large. 



6.4.2 Stress tensor and charge currents 



In appendix IC.21 we have computed the boundary stress tensor corresponding to this black 
hole solution (by foliating the space into S 3 's at infinity, computing the extrinsic curvature 
of these sections, and subtracting the appropriate counterterms) . At leading order in ^- 



m 



8nG, 

m 
8^Gi 



■7 4 (4 7 5 



"7 



4 7 V + - 



8ttG 5 
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8ttG 5 
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8ttG 5 
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8ttG 5 
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a7 u 
a6 7 6 



(95) 



In a similar fashion, the charge currents on S 3 may be computed from J- 1 = —r^g^A^ ^^, 
where the indices fi, v are tangent to the S 3 x time foliations and and the bulk gauge field 



20 The functions h(v) and its derivatives have simple expressions as functions of bulk parameters. 
Comparing with (|64[) we find 



h{v) 
hi{v) 
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16ttT 4 
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A v is given in the equation (2) of 18j. We find 



Jt 


- 7* 

_ J 2 


— J 3 " 


• J 1 


— J 2 


- f 
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— J 2 


- J* 
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— J 2 


- f 



Q 4 
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8wG 5 




q 4 (96) 
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8ttG 5 
g 

8ttG f 



7 4 6. 



Using (j9"%|) . it is evident that the expressions in (|96p are in precise agreement with the 
predictions (lo"Tj) of fluid dynamics. 

6.5 Black holes with independent SO (6) charges and two equal 
rotations 

The most general (five parameter) black hole solutions with the two angular velocities set 
equal can be found in [17| . The thermodynamics of these black holes was computed in 23 



The black hole solutions depend on the parameters 5i, 8 2 , #3, a, m, r + that are related 
by the equation Y(r) = 0. The thermodynamics of these black holes simplify in the limit 

2,TflS 2 

r + > 1 , — = Hi-l fixed. 
Then solving the equation Y = in this limit, one can express m as 



2m 



'I -a 2 



The various thermodynamic quantities in this limi10 (after multiplying integrals by 
can be summarised by 



2tt 

3 

<*= r+VT ^ (-H- ) n - ^ {- 



We believe that [23| has a typo: (3.10) should read $^ = -^^-{siCi + ^aft(ciSjSk — SiCjCk)). Note that 

jd^ + %d$ 2 so that Q a = ^. 



they also use coordinates tp — <$>\ + 02 and </? = </>i — 2 so that 09^, = ^d ( f >1 + §<9</> 2 so that Q a = ^ 
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These expressions match with (1831) if Ki is identified with Hi — 1, demonstrating perfect 
agreement between black hole and fluid dynamical thermodynamics. 

Translating the estimate for the mean free path into the black hole variables, we find 



(an equation that is valid only in the large r + limit). Notice that / m f p is automatically 
small in the large r + limit, explaining why black hole thermodynamics agrees with the 
predictions of the Navier-Stokes equations in this limit. 

Notice that the fluid mechanical expansion parameter Z m f p appears to differ from the 
expansion parameter of black hole thermodynamics used above, l/r+, by a factor of 
1 / VYii Hi- When the three charges of the black hole are in any fixed ratio a : b : c, 
with none of a, b or c either zero or infinity, it may easily be verified that this additional 
factor is bounded between a nonzero number (which depends on a, b, c) and unity. In this 
case the two expansion parameters - / m f p and l/r + - are essentially the same. 

However when one of the black hole charges (say Ri) vanishes H 2 and/or H 3 can for- 
mally take arbitrarily large values. In this extreme limit / m f p appears to differ significantly 
from the bulk expansion parameter l/r + . However large Hi implies large Ki, a limit that 
we have argued above to be thermodynamically singular. Keeping away from the suspi- 
cious large Ki limit, it is always true that Z m f p is essentially identical l/r+, the parameter 
in which we have expanded the formulas of black hole thermodynamics. 

Finally we emphasise that the black hole studied in this subsection include a large class 
of perfectly nonsingular zero temperature or extremal black holes with finite Ki and large 
r + which perfectly reproduce the predictions of extremal fluid mechanics of § §6.21 

In more detail, the thermodynamical quantities of a general solution in this subsection 
is given in terms of X, Y, Z (defined as in (172]) ) as 



S 
E 



a 2 ) 



^XYZ(l 
N 2 r 4 + (3 + a 2 ) 
AXYZil - a 2 ) 



N 2 ria 



L = ± (9£ 

2XYZ(1 - a 2 ) ' 1 



Ri 



N 2 ri II- Xi 



2^XYZ(l-a 2 ) V X, 



/*<i-*)(i 



XYZ 

From these expressions, together with the formula for temperature in ( 1971) it follows that 
the limit X + Y + Z — > 1 (with none of X, Y, Z zero) is extremal (the temperature goes 
to zero) and non-singular (all thermodynamic quantities are finite and well defined). Note 
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that r + is an arbitrary parameter for these extremal black holes. When r + is large the 
fluid dynamical description is valid. The black holes so obtained are exactly dual to the 
extremal fluid configurations described in § §6.2[ 

6.6 Black holes with two equal large R-charges and third R- 
charge small 



Chong et al. [19] have determined a class of black hole solutions with two SO (6) charges 
held equal, while the third charge is varied as a function of these two equal charges. In the 
large radius limit, it turns out that this third charge is negligible compared to the first two, 
so for our purposes these solutions can be thought of as black holes with two equal SO (6) 
charges, with arbitrary rotations and the third SO (6) charge set to zero. The parameters 
of this black hole solution are a, b, m, r+, s, which are related by the equation X(r+) = 0. 
Black hole formulae simplify in the limit 

. , 2ms 2 
r + ^> 1 and k = — = — fixed, 

r\_ 

in units where the inverse AdS radius g = 1, which leads to 

2m = 7^(1 + A;) 2 . 

Multiplying all thermodynamic integrals in [l^] by = 2j ^-, in this limit, the ther- 
modynamics can be summarised by 

Qi = a , Q2 = b , T = — , 

7T 



Cl = C2 = 7r 2Vfc, c 3 ~o(J-) 

_ 7r 2 jV 2 y 4 T 3 (l + k) 2 
gc ~ 8(l-a 2 )(l-6 2 ) • 



(99) 



Note that (3 and R3 are subleading in r + . These formulae are in perfect agreement with 
( 1861) if we identify 

k — k . 

From the expression for the temperature, it follows that all extremal or zero temperature 
black holes have r + = 0. Consequently all extremal black holes (of the class of black holes 
described in this subsection) are singular, dual to the fact that the fluid mechanics has no 
thermodynamically nonsingular zero temperature solutions. 

Translating the estimate for the fluid dynamical mean free path into the black hole 
variables we find (assuming r + 3> 0) 

/mfp ~ 3r + (l + K ) • 

It follows that the fluid dynamical expansion parameter is essentially the same as l/r + , 
provided we stay away from the thermodynamically suspect parameter regime of large k . 
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6.7 Black holes with two R-charges zero 

The solution for the most general black hole with two R-charges set to zero relevant solution 
has was presented in 2(|. The parameters of this black hole are xo,m, 6, a,b related by 
X(x ) = 0. 

The thermodynamics of these black holes simplifies in the limit 

x > 1 , y = ^ s /x^5 fixed, 
in units where g — 1, which leads to 



2m 

This gives an upper bound on y: y < 1. 

Multiplying all thermodynamic inte; 
modynamic formulae can be summarised by 



i-y 



2\ 



Multiplying all thermodynamic integrals in 2(| by = in this limit, the ther- 



Qi — a , VL 2 = b . 



y/^(2 - y 2 ) 2ixTy^Y^y~ 2 
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InZ 



2nyn—yZ v 2-y 2 (100) 

xl /2 nN 2 47r 4 AT 2 T 3 (l - y 2 ) 



SC 2^1—^(2 - y 2 )(l - a 2 )(l - b 2 ) (1 - 0?)(1 - ^)(2 - y 2 Y 

2 

Upon identifying k = jz^z, we find perfect agreement with (|89|) . Under this identification, 
the expression for temperature becomes 

rp v^o( 2 + K ) 



27TA/1 + K 

As in the subsection above, it follows immediately from this equation that the black hole 
temperature vanishes only for the singular black holes with xq = 0. This matches with the 
fact that there are no nonsingular extremal fluid dynamical solutions in this case. 

The fluid dynamical mean free path may be evaluated as a function of bulk parameters 

as 

^mfp ~ = ■ 



3Va;o(l + k) 



Note that Z mfp is small whenever ^/xo = r + is large, an observation that explains the 
agreement of black hole thermodynamics in the large r + limit with the Navier-Stokes 
equations. In more generality we see that Z m f p is essentially the same as l/r + , provided we 
keep away from the thermodynamically suspicious parameter regime of k large. 
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6.8 Extremality and the attractor mechanism 

As discussed in the previous subsections, there exists a duality between extremal large 
rotating AdS black holes on one hand and the extremal configurations of the fluid dynam- 
ics on the other. This implies that the thermodynamic properties of these large rotating 
extremal black holes are completely determined by the corresponding properties of large 
static extremal black holes. As an application of this observation, let us recall the sugges- 
tion 3lM33l| that the attractor mechanism for black holes implies the non-renormalisation 
of the entropy of all extremal configurations, as a function of the 't Hooft coupling A. It 
follows immediately from the fluid mechanical description at large charges, that were any 
such non-renormalisation theorem be proved for static extremal configurations, it would 
immediately imply a similar result for rotating extremal configurations. 

6.9 BPS bound and supersymmetric black holes 

All solutions of IIB supergravity on AdS$ x S 5 , and all configurations of M = 4 Yang-Mills 
on S 3 obey the BPS bound 

E > L x + L 2 + ^Ri = Li + L 2 + 3R. (101) 

i 

Within the validity of the fluid dynamical approximation, described in this paper, 

E - L x - L 2 = 2tt 2 T 4 A 3 / + ^ +^-^2 (102) 

(l+wi)(l + w 2 ) 1 ; 

notice that the RHS of this equation is positive definite. The BPS bound is obeyed provided 

3 + ^+^-a;^ 
(l + wiXl + Wa) " K 1 

Plugging in the explicit expressions for A and Oj from (1771) . we find this condition is satisfied 
provided 

r+ = > Vgi±Ml±^ . (104 ) 

2 - K (1 + K){3 + UJi + UJ 2 - UJ\UJ 2 ) 

The RHS of (I104p is of order unity. It follows that (11041) is saturated only when r + of unit 
order. It follows that when r + ^> 1 (so that fluid dynamics is a valid approximation) the 
BPS bound is always obeyed as a strict inequality. Supersymmetric black holes are never 
reliably described within fluid mechanics@ The extremal black holes with large horizon 



22 Although it is possible to make the energy of supersymmetric black holes parametrically larger than 
their entropy, this is achieved by scaling either uj\ or u>2 to unity with r + kept at unit order. It is easy to 
verify that in this limit the local, rest frame mean free path of the fluid is of unit order in regions of the 
S 3 and so fluid mechanics may not be used to describe these configurations. 
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radius, that are well described by fluid mechanic^ (see the previous subsection) are always 
far from supersymmetry. 

We have noted above that a large class of extremal configurations in strongly interacting 
Yang-Mills - all those that admit a fluid dynamic description - are not BPS. This is in 
sharp contrast with the results of computations in free Yang-Mills theory, in which all 
extremal configurations are supersymmetric [34|. This difference is related to the fact, 
noted previously, the divergent mean free path prevents a fluid mechanical description 
from applying to free theories. A practical manifestation of this fact is that the function 
h{y), which appears in the analysis of free Yang-Mills in equation (5.2) of js^], and plays 
the role of r + in our discussion here, is always of order unity for all allowed values of the 
chemical potential, and so can never become large. 

6.10 Fluid dynamics versus black hole physics at next to leading 
order 

As we have explained above, the formulae for all thermodynamic charges and potentials of 
black holes of temperature T and chemical potentials Ui, in AdS§ x S 5 , may be expanded 
as a Taylor series in l/r+ ~ l m tp(T, z/j). As we have verified above, for every known family 
of large AdS black holes, the leading order results in this expansion perfectly match the 
predictions of the Navier-Stokes equations. Higher order terms in this expansion represent 
corrections to Navier-Stokes equations. In this subsection we investigate the structure of 
these corrections. 

Let us first investigate the case of black holes with at least one SO (6) charge set equal 
to zero (the black holes studied in § §6.61 and § §6.7p . It is not difficult to verify that the first 
deviations from the large radius thermodynamics of these black holes occur at 0(1/7+) ~ 
^mfp- This result is in perfect accord with naive expectations from fluid mechanics. As we 
have explained above, the fluid dynamical configurations presented in this paper are exact 
solutions to the equations of fluid mechanics with all one derivative terms, i.e. to the first 
order in / mfp . In general we would expect our solutions (and their thermodynamics) to be 
modified at 0(1^ ), exactly as we find from the black hole formulae. 

However when we turn our attention to black holes with all three SO (6) charges nonzero 
we run into a bit of a surprise. It appears that the thermodynamics (and stress tensor and 
charge currents) of these black holes receives corrections at order 0(l/r + ) ~ / m f p . This 
result is a surprise because, for the reason we have explained in the previous paragraph, 
we would have expected the first corrections to our fluid mechanical configuration to occur 

ate>(£ fp ). 

We do not have a satisfactory resolution to this puzzle. In this subsection we will 
simply present the expressions for the first order corrections to black hole thermodynamics 
in a particular case (the case of black holes with all SO (6) charges equal), and leave the 

23 Note that the 'physical' radius (Area) 1 ' 3 of the black hole is distinct from the parameter r + which 
determines the validity of fluid dynamics. The physical radius can be made arbitrarily large, nevertheless 
fluid mechanics is only valid if r + is large. 
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explanation of these formulae to future work. 

As we have mentioned above, the thermodynamics of a charged rotating black hole in 
AdSs x S 5 with three equal charges and two different angular momenta can be found in 
181 ] . To calculate next to leading order (NLO) corrections to the thermodynamics of large 
black holes, we systematically expand the thermodynamic quantities. 

We find it convenient to shift to a new parametrisation in which there are no NLO 
corrections to the intensive quantities. This allows us to cast the NLO corrections entirely 
in terms of the intensive quantities. The parameters we choose are related to the parameters 
in 



181 ] in the following way 
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(105) 



r + = £ + ^Kuj a uj b , 
q = + 3n£ 2 uj a uj b . 

In terms of these parameters, the intensive quantities can be written as 

n a = u a + o ' 1 



Mb 
T 



oo b + 

2-K 

2tt-\/k 



+ 
-O 



(106) 



where we have calculated up to NLO and confirmed that the intensive quantities do not 
get corrected in this order. 

This in turn means that the new parameters can be directly interpreted in terms of the 
intensive quantities. 



na + o[U, uj b = n b + o[i 2 mip ] 



where Z mfp ~ ^p. 
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Now, we calculate NLO corrections to the extensive quantities in terms of the new 
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parameters. 
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where G 5 = irR\ dS /(2N 2 ) is the Newton's constant in AdS 5 



In particular, the subleading terms can be isolated and written as 
AS = 0, 
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7 Comparison with black holes in AdS^ x S 7 and AdSV x 

In this section we compare solutions of rotating fluids of the M5 or M2 brane conformal 
field theory on S* 2 or S 5 to the classical physics of black holes in M theory on AdS^ x S 7 
and AdS-j x S 5 respectively. Our results turn out to be qualitatively similar to those of 
the previous section with one difference: the puzzle regarding the next to leading order 
agreement between fluid dynamics and black hole physics seems to be absent in this case. 
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7.1 Predictions from fluid mechanics 



The equations of state of the strongly coupled M2 and M5 brane fluids were computed 
from spinning brane solutions in 35[. Our parameters are related to theirs by = lf/r 2 H . 



7.1.1 M2 branes 



We define our R-charges to be half of the angular momenta of 35j to agree with gauged 
supergravity conventions. The equation of state is 
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where k — 1 . . . 4. 

The stress tensor and currents are given by (1501) and (15T1) with 
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3(3 + 2 £\ + E j<fc KjKk - Uj Kjf ' 

4 7 r 2 (2iV) 3 / 2 n j (l + ^) 3/2 
3(3 + 2 V. Kj + V < fc « iKfc - [L Kj) 2 ' 



vr(2iV) 3 / 2 n i (l + ^) 3/2 



3(3 + 2 Y,j *j + T,j<k K i K k ~ Eli «j) 



The thermodynamics can be summarised by 

47rrn J (i+« i ; 



InZ. 



(3+2 £\ Kj - + ^ i<fc « iKfc - rij v 1 + «• 

16vr 3 (2iV) 3 / 2 r 2 n 7 (l + s) 5/2 



3(1 - q 2 )(3 + 2 £ . «i + E s<k ^ - n- ' 

The mean free path in fluid dynamics is given by 



'mfp 



3 



B 



Q=0 (d - l)AnTA 

V— 

4^1 + 



87rTlI,-(l + «i) 



1 

8^T 



(109) 
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This simplifies when the charges are pairwise equal, K3 = K\ and K4 = K2. In this case, 
with i = 1,2: 



InZ 



4ttTIL(1 + «, 



(3 + J2j i; j - \ [ j K j) V 1 • 

167r 3 (2A^) 3 / 2 T 2 n ? (l + s)' 



3(1- n*)(3 + n,*;) 3 ' 



113) 



and the mean free path becomes 

+ j_ (114) 

87^11,(1 + SiT^l + Ki J' 

It is evident that the thermodynamic equations of state listed above allow a set of 
extremal fluid configurations very similar to those discussed in § §6.21 The analysis of § §6.21 
can be easily extended to fluids on S 2 . 



7.1.2 M5 branes 



We define our R-charges to be twice the angular momenta of [35j to agree with gauged 
supergravity conventions. The equation of state is 

64vr 3 Ar 3 n ? (l + S-) 4 



h(u) 



27^(1 + 



hi(v) 



(3+E^i-n,^) vi + ^ 

128vr 2 iV 3 n,(l + s) : 

3(3 • v ; k ; ii^,-;; 



where i = 1,2. 



The stress tensor and currents are given by (I50I) and (|5T|) with 

647r 3 iV 3 n i (i + '%) 4 

128vr 3 iV 3 Jl^l + Kjf 
3(3 • ll^/o) 5 

i287r 2 iv 3 ni(i + ^) 3 

3(3 • V ; /,- ; I I y / , 5 ' 
The thermodynamics can be summarised by 

47TTiyi + Ki) 



.4 



5 



InZ. 



(3+2 £\ Ki + ^ i<fc KjK k - n 3 - «i) v 1 + ^ 

64vr 6 iV 3 T 5 n^i + S) 4 

sriaa-^Ks+Ei^-n,-^) 3 ' 



(115) 



(116) 



(117) 
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The mean free path in fluid dynamics is given by 



'mfp 



s 



AtxE 



B 



(s + y , Kj - Uj '<•„ 



(d - i)AttTA iottT n (1 + 



IOttT 



^ 1 + K 



- 1 



In the case that the three rotation parameters are equal, Qi = Q2 = ^3 
7 = (1 - VI 2 )- 1 / 2 and 



;ii8) 



Q, we have 



47rTn(i + «i 



(3 + 2 V . Kj + V . <fc ^-/Cfc - Uj K j) V 1 + K 



InZ. 



64yr 6 iV 3 T 5 n 7 + 



;ii9) 



3 (l_ n2) 3 ( 3 + E n .) 



3 ' 



It is evident that the thermodynamic equations of state listed above allow a set of extremal 
fluid configurations very similar to those discussed in § §6.21 The analysis of § §6.21 can be 
easily extended to fluids on S 5 . 



7.2 Black holes in AdS4 with pairwise equal charges 



The relevant solution was found in [22j. Its thermodynamics have been computed in [23 
We consider the limit of large r + with = hi fixed. In this limit m can be written as 

3 

m = r -±(l + kl ) 2 (l + k 2 ) 2 , 

and therefore s,- ~ — . 

After multiplying integrals by -£f^- = ^ 2N ^ 1 , the thermodynamic quantities can be 
expressed as 

T 'M3 • V /Q ■ - 1 1 f /Q) 

T= ^ , it — a, 

4:71 

Ci = C3=4 " t (3Te^TW' C2 ^ C4 (3+E,%-n,^)' ( 12 °) 

167r 3 (2iV) 3 / 2 T 2 / n,(l + ^) 2 \ 1 



InZ, 



I ("? _i_ . _ FT I- AS I 1_ fl 2 



(3+ea-il%: 



If one identifies Ki = hi, then these formulae match with ( 11131) . It is not difficult to verify 
that the first corrections to the thermodynamical equations above occur at 0(l/r\). 

It is clear from ( 11201) that the black holes of this subsection admit a zero tempera- 
ture (extremal) limit with nonsingular thermodynamics at any every value of r + . These 



43 



extremal black holes are dual to extremal solutions of fluid dynamics analogous to those 
described in the previous section in the context of TV = 4 Yang-Mills. 

The fluid dynamical mean free path may easily be computed as a function of black hole 
parameters. From (1114p we find 



/mfp ~2r + n,(l + ^)- 

As in the previous section, the Z m f p ~ l/r + away from thermodynamically suspect limits 
of parameters. 



7.3 Black holes in AdSy with equal rotation parameters 

The relevant solution was found in [2lj]. Its thermodynamics have been computed in j23|@ 
We set the parameter g in 23] to be unity and consider the limit 

p + ^> 1 , and Hi = 1 H tA- fixed, 

P+ 

where i=l,2. In this limit, the parameter m is given by 

In this limit, after multiplying integrals by = the thermodynamics can be 

summarised by 

Q = a , T 



^ { IL vh, 



Cl = ^>:A ii.// ■ 6 = 2Tr 2 £j X-n^ • (12i; 



InZ 



647r 6 iV 3 T 5 / UjHf 



sc 3(1-^)3 ^(2E J ^-n i ^i) 6 



These formulae agree with (11191) upon identifying «j = H^ — l. The first corrections to these 
thermodynamical formulae occur at 0{l/rV). Using this identification we can rewrite the 
expression for the temperature as 



2tt I ^ 1 + « 



24 We believe that [23[ has the following typos: equation (4.7) should read 
7r 3 (r 2 +a 2 )JJl Y' 2m Si , 
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It follows that the black holes studied in this subsection admit smooth extremal limits at 
any value of p + . Extremal black holes with large p + (and with no Ki arbitrarily large) are 
dual to extremal solutions of fluid mechanics. 

Expressing the fluid mechanical mean free path d 1 1 8 f) as a function of black hole pa- 
rameters we find ^ 

ZmfP 5p+ Uj V 1 + K i ' 
Once again Z m f p ~ away from thermodynamically suspect limits. 



8 Discussion 

As we have explained in this paper, the classical properties of large black holes in AdS 
spaces enjoy a large degree of universality, summarised by (j2J). However the reasoning that 
led to flSJ) applies equally to all classical theories of gravity, not just to those theories that 
are governed by the two derivative effective action. For instance, Af = 4 Yang-Mills theory 
at finite A is dual to IIB theory on AdS§ x S 5 of finite radius in string units. Even though 
thermodynamics of black holes in this background will receive contributions from each of 
the infinite sequence of a' corrections to the Einstein action, we expect ([2]) to be exact in 
the large horizon radius limit @ 

We find it particularly interesting that (at least in several particular contexts) our fluid 
dynamical picture applies not just to non-extremal black holes but also to large radius 
extremal black holes. This fact might allow us to make connections between our approach 
and the interesting recent investigations of the properties of extremal black holes. In 



particular, Astefanesei, Goldstein, Jena, Sen and Trivedi [33( have recently argued that 
the attractor mechanism applies to rotating extremal black holes, and have derived a 
differential equation that determines the attractor geometry (and gauge field distribution, 
etc.) of the near horizon region of such black holes. It would be very interesting to 
investigate the connection, if any, between these rotating attractor equations and our 
equations of rotating fluid dynamics. 

It would be conceptually simple (though perhaps technically intricate) to compute the 
spectrum of small fluctuations about the fluid dynamical solutions presented in this paper. 
This spectrum should match the spectrum of the (lowest) quasinormal modes about the 
relevant black holes (the decay of fluid fluctuations due to viscosity maps to the decay of 
quasi normal modes as they fall into the black hole horizon). It would be interesting to 
check if this is indeed the case. 

It would be interesting to better understand, purely in bulk terms, why our proposal 
works. Roughly speaking, it should be possible to understand the metric of a black hole 
in global AdS and in the large radius limit, as a superposition of patches of the metric of 



25 Away from the supergravity limit, the mean free path Z m f p — v j p is expected to be given by f(X)s/p 
where /(A) is a monotonically decreasing function that interpolates between infinity at A = to unity at 
infinite A. Thus the condition for the validity of fluid mechanics is modified at finite A; in the uncharged 
case, for instance, it is T ^ /(A). 
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black branes of various different temperatures and moving at various different velocities, 
where the temperatures and velocities are given by the solutions to the fluid dynamical 
equations presented in this paper. It would be interesting and useful if these words could 
be converted into the first term of a systematic approximation procedure to generate black 
hole solutions in AdS spaces in a power series in l/r+. Such a construction would constitute 
a bulk derivation of the boundary Navier-Stokes equations (and corrections thereof). 

Relatedly, it would be interesting to ask if there are any gravitational interpretations of 
the local properties of fluids in our solutions. For instance, fluid mechanics yields a sharp 
prediction for the velocity and entropy density of the fluid as a function of position on the 
sphere. We have not yet been able to verify these predictions, because we do not know 
what gravitational construction we should compare them to. The entropy of the fluid is an 
integral over the boundary. The entropy of the black hole is an integral over the horizon. 
Perhaps there exists a natural map from the horizon to the boundary that allows one to 
convert horizon densities to boundary densities and vice versa. Such a map (for which we 
have no conjecture) would permit a gravitational prediction s, the local entropy density of 
the fluid. 

The fluid velocity is another quantity for which it would be useful to have a gravitational 
definition. We do not really have a serious proposal for such a definition: nonetheless, in 
the next few paragraphs we outline a caricature proposal, in order to give the reader a 
sense of the types of relations that might exist (we emphasise that we do not have any 
physical reason to believe that this caricature has any truth to it). 

In the black hole solutions, there is one special Killing vector, K — d t + Vt a d^ a) that is 
also the null generator of the horizon. It has the norm 



\K\\ 2 = -K^K^ 



r 2 7 2 at the boundary, 
at the horizon. 



If we were to normalise it with respect to the metric of the conformal boundary, the result 
would be K —jK. This could be identified with the fluid velocity m m . However, as 7 is 
not constant, K is not a Killing vector. It also seems unnatural to use a normalisation 
factor that depends on 9 but not r. Nonetheless, this much maligned vector field has an 
interesting property. 

Recall that black hole temperature and chemical potentials can be computed from the 
formulae 



2tt 2tt 



horizon 



horizon 



If one were to replace K with K in the formulae above, one would obtain T and //j, the 
local temperature and chemical potentials of the fluid. 

We end this paper by reminding the reader that, while our proposal has passed many 
checks, our work has left one significant puzzle unresolved. While the thermodynamics 
and stress tensors of uncharged rotating black holes in every dimension, plus all known 
black holes in AdS-j x S* 4 and AdS± x S 7 , deviate from the predictions of the Navier-Stokes 
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equations only at second order in / m f p , the situation is more complex for black holes in 
AdSz, x S 5 . In this case, black holes with at least one SO (6) charge equal to zero also agree 
with the results of the Navier-Stokes equations up to C?(Z^ fp ). However the thermodynamics 
of rotating black holes with all SO (6) charges nonzero, appears to deviate from our fluid 
mechanical predictions at C(/ m fp) (see § §6. 10j) . 

We consider this a significant puzzle as our fluid dynamical configurations solve the 
Navier-Stokes equations including 0(/ m f p ) dissipative contributions. Moreover, in appendix 
|A]we have checked by direct enumeration that all possible vectors and traceless symmetric 
tensors that transform homogeneously under conformal transformation and contain a single 
derivative simply vanish on our solution, so it is difficult to see how any one derivative 
modification to the equations of fluid dynamics could help resolve this puzzle 1^1 Once this 
puzzle is resolved it would be interesting to attempt to reproduce the C(/^ fp ) corrections to 
black hole thermodynamics from appropriate additions to the equations of fluid dynamics. 
It is perhaps worth emphasising that black holes in AdS represent exact (to all orders in 
^mfp) solutions to a dynamical flow. A detailed study of these solutions might lead to new 
insights into the nature of the fluid dynamical approximations of the high energy regime 
of quantum field theories. 
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Appendices 

A Conformal fluid mechanics 

Consider a conformal fluid in d dimensions. We seek the conformal transformations of 
various observables of such a fluid. To this end, consider a conformal transformation which 

26 It has been suggested that certain pathologies in relativistic fluid dynamics lead to the breakdown of 
the derivative expansion [36| (See also Q and the references therein) . As any such pathology should apply 
equally to two charge and three charge black holes, we find it difficult to see how this issue could have 
bearing on our puzzle. We thank S. Gupta and H. Liu for discussions on this issue. Another possibility 
is that the formulae of black hole thermodynamics receive corrections - perhaps from Wess-Zumino type 
terms - that are nonzero only in an even dimensional bulk (and so in ten but not in eleven dimensions) 
and only when all charges are nonzero. We thank O. Aharony for suggesting this possibility. 
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replaces the old metric g^ v with g^ given by 



_ 20- 

9 fxv ^ 9fj.ui 



The Christoffel symbols transform as 



pi/ 



Let m m be the four- velocity describing the fluid motion. Using g^u^u" = g iiv u il u v = - 
we get m m = e~^u^. It follows that the projection tensor transforms as P^ u = g^ u + u^u v 
Q- 2( t>P^ u . The transformation of the covariant derivative of uF" is given by 



122) 



The above equation can be used to derive the transformation of various related quan- 
tities 



& = VX = e" 



vPVpU" = e 



-2<f> 



d-1 



-$pV = e -30^ 



(123) 



Further, the transformation of the temperature and the chemical potential can be written 
as T = e~^T and fi = e~^jl. The transformation of spatial gradient of temperature 
(appearing in the Fourier law of heat conduction) is 



P^(d v T + a u T) = e- 3 *P» u (d v T + a v T) ■ 

The viscosity, conductivity etc. scale as k = e~^ d ~ 2 ^k , r\ = e~^ d ~ 1 '^'fj 
Dij = e-W+Dy. 

For a fluid with c charges, there are 2c + 2 vector quantities involving no more than a 
single derivative which transform homogeneously^!. They are 



fj, ^ — e and 



d^, d u T + ( a 



x) 



d- 1 



Un T, u^vfdvVi and ( u a d a T + 



■d 



d-1 



In the kind of solutions we consider in this paper, all of them vanish except w M . 

The transformation of the stress tensor is T^ v = Q- ( , d + 2 )<t> r p^ i/ ) f r0 m which it follows that 



27 In the following analysis, we will neglect pseudo-tensors which can be formed out of e M „.... Additional 
tensors appear if such pseudo-tensors are included in the analysis. 
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So, for the stress tensor to be conserved in both the metrics, it is necessary that T^ v is 
traceless. 

To consider the possible terms that can appear in the stress tensor, we should look at 
the traceless symmetric second rank tensors which transform homogeneously. The tensors 
formed out of single derivatives which satisfy the above criterion are easily enumerated. 
For a fluid with c charges, there are 2c + 4 such tensors and they are 

uV + ^g^, a^, q»u v + q v u^ (u a d a T + j^j?) (w*u u + ^ 



(124) 

1 a^ v ( 1 \ 
- (u^u, + u x d^i) - *— u a d a Vi and u a d a ^ u»u v + -g^ . 

2 ' d \ a J 



Among these possibilities, the stress tensor we employ just contains the tensors in the first 
line. It can be shown that the other tensors which appear in the above list can be removed 
by a redefinition of the temperature etc. Even if they were to appear in the stress tensor, 
for the purposes of this paper, it suffices to notice that all such tensors except u^u" + \g^ v 
vanish on our solutions. Hence, they would not contribute to any of the thermodynamic 
integrals evaluated on our solutions. 



B Free thermodynamics on spheres 

In ( |55l) above, we have presented a general expression for the grand canonical partition 
function for any conformal fluid on a sphere. In this appendix , we compare this expression 
with the conformal thermodynamics of a free complex scalar field on a sphere. 

Strictly speaking, the fluid dynamical description never applies to free theories on a 
compact manifold, as the constituents of a free gas have a divergent mean free path (they 
never collide). Nonetheless, as we demonstrate in this subsection, free thermodynamics 
already displays some of the features of fl55|) - in its dependence on angular velocities, for 
example - together with certain pathologies unique to free theories. 

Consider a free complex scalar field on S d ~ 1 x time. This system has a U(l) symmetry, 
under which has unit charge and <fi* has charge minus one. We define the 'letter partition 



function' [37|] Z\ et as Tr exp [— (3H + vR + /3Q a L a ] evaluated over all spherical harmonic 



modes of the scalar field 



,d-2 



1 — e~ 2/3 

>iet = (e" + e">-^ ( v n Li (i- e -^)(l-e-^) > ll2f,) 



(this formula, and some of the others in this section, are valid only for even d; the gen- 
eralisation to odd d is simple). We will now examine the high temperature limit of the 
grand-canonical partition function separately for v = and v ^ 0. 
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B.l Zero chemical potential: {y = 0) case 

The second quantised partition function, Z gc for the scalar field on the sphere is given by 



exp Ij2 

\ N 



Z gc = exp [ ^ ) . (126) 



For small /?, we have 



J lct 



/3 d - 1 IL(i - «2) 

It follows thafi 



Upon identifying Vdh\ u=0 = 4£(<i), we find that (11261) is in perfect agreement with ([55 



B.2 Nonzero chemical potential: {y ^ 0) case 

The high temperature limit of the thermodynamics of a free, charged, massless field is 
complicated by the occurrence of Bose condensation. This phenomenon occurs already 
when uj a = 0; this is the case we first focus on. 

It is useful to rewrite the letter partition function as 

Z let = (2 cosh u) e"^ ™(N)e~ m } (128) 

N 

where m(N) w 2N d ~ 2 /(d — 2)! for N ^> 1. The logarithm of the grand canonical partition 
function may then be written as a sum over Bose factors (one per 'letter') 

In Z gc = - m ( N ) M 1 - e^ 7 ^- 2 ^) + ln(l - e -WHd-2)/2)-vfl . (129 ) 

N 

The total charge in this ensemble is given by 

R= -^ lnZ *> = T, m(iV) ( e/W(d-2)/2)-» - 1 " e ^+(^)/2) + ._i ) ' ( 13 °) 

In order to compare with fluid dynamics, we should take (3 to zero while simultaneously 
scaling to large R as R = ^Jtt with q held fixed. As we will see below, in order to make the 
total charge R large, we will have to choose the chemical potential to be large. However 
it is clear from (11291) that \u\ < (3(d — 2)/2. Consequently, the best we can do is to set 

28 This formula has been derived before in many contexts, for example [38} have derived this in d = 4 
and compared it with the thermodynamics of black holes in AdSs . 
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v = f3((d — 2)/2) — e where e will be taken to be small. We are interested in the limit when 
3 is also small. We may approximate (j!30p by 



1 A / 1 1 



fid-! e e P(d-2)-e _ 1 + (^JV+e _ I e/ a(JV+(d-2))-e - 1 J ' ^ 131 ^ 

The only solution to f 1 1 3 1 j) is 

3 d ' 1 
1 

Substituting this solution into the partition function, we find 

\nZ q = ^-(l + 0(3)). (132) 

Consequently, to leading order the partition function is independent of the charge q ! What 
is going on here is that almost all of the charge of the system resides in a Bose condensate 
of the zero mode of the field <fi. This zero mode contributes very little entropy or energy to 
the system at leading order in /3@ At high temperatures, the zero mode is simply a sink 
that absorbs the system charge, leaving the other thermodynamic parameters unaffected. 

Upon generalising our analysis to include angular velocities, we once again find that the 
leading order partition function (in the limit of high temperatures and a charge R = q/fi^ 1 ) 
is independent of q and in fact is given by ( 1127R . Consequently, there is a slightly trivial (or 
pathological) sense in which the thermodynamics of a free charged scalar field agrees with 
the predictions of fluid mechanics - we find agreement upon setting h{y) to a constant. 



C Stress tensors from black holes 

According to the usual AdS/CFT dictionary, the boundary stress tensor on S d ~ l , corre- 
sponding to any finite energy solution about an AdS^ background of gravity may be read 
off from the metric near the boundary, using the following procedure (3-lo|. First we foliate 



the spacetime near the boundary into a one parameter set of d geometries, each of which 
is metrically conformal to S d ~ l x R, to leading order in deviations from the boundary. 
We will find it convenient to use coordinates such that the leading order metric in the 
neighbourhood of the boundary takes the form 

ds 2 = -r 2 dt 2 + ^dQ 2 ,^ + ^- (133) 

(here r = oo is the boundary; this metric has corrections at sub leading orders in -?). 

In these coordinates, our foliation surfaces are simply given by r = const. We next 
compute the extrinsic curvature 0^ = — V ll n u on these surfaces, where n u is the unit 



29 In particular, the contribution of the zero mode to the energy is proportional to the charge, which is 
suppressed by a factor of [3 relative to the contribution to the energy from nonzero modes. 
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outward normal to these surfaces. The boundary stress tensor for the dual field theory on 



a unit sphere is given by [3 10 



r D-l 



W u = lim — - (0* - ^0) , (134) 

where the coordinates /i, v go over time and the angles on S . 

The stress tensor as defined above will contain some terms which are independent of 
mass and charge of the black hole. These are the terms that are nonzero even on the 
vacuum AdS background and they diverge in the limit r — > 00. These terms are all 
precisely cancelled, up to a zero point Casimir energy, by counter terms presented in §2 of 
(3j. We will simply ignore all such terms below; consequently, the stress tensors we present 
in this paper should be thought of as the field theory stress tensors with the contribution 
from the Casimir energy subtracted out] 3 "*! 

In order to compute the stress tensor in (11341) . we must retain subleading corrections 
to the metric in (1133p . However, only those corrections that are subleading at (9(-d=t) 
(relative to the leading order metric in (11331) ) contribute to (j!34p . 

In order to compute the stress tensor corresponding to two classes of black hole solutions 
below, we adopt the following procedure. First, we find a coordinate change that casts the 
metric at infinity in the form (11331) at leading order. Next we compute all the subleading 
corrections to the metric at order 0( Finally, we use these corrections to compute 

the extrinsic curvature 9^ and then IT^ using (11341) . 



C.l Stress tensor from rotating uncharged black holes in AdS^j 
C.l.l D = 2n + 1 

The most general rotating black hole in an odd dimensional AdS space is given by equation 
(E.3) of |l5| (we specialise to odd dimensions by setting the parameter e in that equation 



to zero). The metric presented in [15j uses as coordinates 



1. The n Killing azimuthal angles fa along which the black hole rotates. These may be 
identified with the coordinates fa in §0] of our paper. 

2. n other unspecified variables (called 'direction cosines') /ij subject to the constraint 
Y^i^i = 1- These may be thought of as the remaining n — 1 coordinates on S* 2 ™ -1 . 

3. The radial variable r and timelike variable t. 



30 We rescale the stress tensor of [3] by a factor of 87r g in order that the energy of our solutions is given 
by J */7 IIq with no extra normalisation factor. 

31 Recall that the full stress tensor of a general d dimensional conformal field theory is not traceless on 
an arbitrary manifold; however the trace is given by a function of the manifold curvature independent of 
the field theory configuration. It follows that our stress tensor with Casimir contribution subtracted must 
be traceless, as indeed it will turn out to be. 
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In order to cast the metric of 15j into the form (I133P near the boundary, we perform 
the following change of coordinates 

r =2_^ — , r^ i (l-o i ) = (r +aj/i i . (135) 

i=i a * 

Note that 

n ra 
i=l i=l 

This equation may be solved by writing fa as functions of the n — 1 variables ^ (which 
may then be identified with the coordinates used in §1]) 

/i-i \ 
/ij = | cos 2 9j J sin 2 ^ . 

In these coordinates, the metric in the neighbourhood of r — ► oo becomes 

df 2 



ds 2 



1 + r 2 )dt 2 + ^L- + r~ 2 ^(d/i 2 + /2 2 d0 2 



i=i 



+ J^7 2( " +1) dt 2 + J^7 2 "df 2 (136) 



fln-2 1 f2n+2 



t=i »=ij=i 

where we have retained all terms that are sub leading up to 0(-pj=r) compared to the metric 
of pure AdS. Here 7 ~ 2 = 1 - £?=i a*tf and Zti ftdfi = (iK cos 2 ^) d£ 2 as in 

Note that this metric separates into two parts; the first piece (on the first line of (j!36p ) 
is the metric of pure AdS space while the terms of the remaining lines represent correction 
proportional to the mass m. 

The normal vector is given by n f = 7= (with all other components zero). As our 

metric contains no terms that mix r with other coordinates at leading order, this is the 
same as 



f 1 ~L 1 
n = — — = r 1 + — 

V9ff V r 



m 2n 



f2n^ 



Since the normal vector has only the r component and since the component is a function 
of f only, to compute the extrinsic curvature tensor Q v x one needs only those components 
of T that are of the form Y v x ~. The Christoffel symbols (that are relevant for the calculation 
of the stress-tensor) as calculated from this metric up to the first subleading term in f are 
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given by 



f 2 + 1 



f / 2nm 2 (n+i)\ r & _ 2nmai 2(n+1) 



t __2nmo^ 2(n+1) , 1 / 2nmaiajfij 2(n+ i)\ (137) 

L ^f- ~ 2n +l I 1 fyr — ~ ~2n ' J V ' 



T e _1 

r 

The extrinsic curvature Q v x in this case is given by 

©A = "^X + 0- ' 



f2n+l 



Ignoring all terms in G^ that are independent of mass (for the reasons explained in the 
introduction to this appendix) we find 

e ^ = ~d=i a ^ °Z = ~5=i a ^ 3 b^j) ( 138 ) 

, _ (D - l) mi D+1 9i _ m 7 ( D ' 1 ) 

Here the n has been replaced by ^y^- It may easily be verified that Q\ is traceless and 
therefore the stress tensor is also traceless according to the definition (I134p . Raising one 
index in G by asymptotic AdS metric, normalising it appropriately and then taking the 
large f limit one can derive the stress tensor as given in (|(d51) . 

C.1.2 D = 2n + 2 

The computation of the boundary stress tensor for the most general uncharged rotating 
black hole in even dimensional AdS spaces is almost identical to the he analysis presented 



in the previous subsection. Once again the metric is given in equation (E-3) of [15j , where 
we must set e to 1 to specialise to even dimensions. The coordinates of the black hole 
solution are similar to those described in the previous subsection, except that there are 
ri+1 coordinates fa restricted by a single equation J2 a A*! = 1- Repeating the computations 
described in the previous subsection, our final result is once again simply (|138|) . In summary 
(11381) is correct no matter whether D is odd or even. 

C.2 Black holes in AdSs with all R-charges equal 

In this subappendix, we compute the boundary stress tensor for a class of charged black 
holes, namely for black holes in AdSs with all three R-charges equal. Our computation 
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will verify the striking prediction of §3] that the functional form of this stress tensor is 
independent of the black hole charge in the fluid dynamical limit. 

The metric for rotating black holes with all R-charges equal is given by (equation (1) 

of (3) 



ds 2 = A- e [{l + y 2 )p 2 dt + 2qu}dt + 2quu + f_ ( Aedt _Y + P 2 dy 2 



S a S 6 p 2 p 2 p 4 V S aSfe J Ay 

p 2 d6 2 y 2 + a 2 . 2 ~ 2 y 2 + b 2 2 ~ 2 
H — 1 — sin 0d0 H — — cos Vdip , 



(139) 



where 



(140) 



A _ {y 2 + a 2 ){y 2 + b 2 )(l + y 2 ) + q 2 + 2abq ^ 

y y 2 
p 2 = y 2 + a 2 cos 2 (9 + b 2 sin 2 , 
A 9 - =1 - a 2 cos 2 - 6 2 sin 2 , 
S a =1 - a 2 , 
£ fe =l-o 2 , 

/ =2mp 2 — q 2 + 2abqp 2 , 

v =b sin 2 0d0 + a cos 2 0dV> , 

uj =asm 0— — ho cos 9-—. 

This metric takes the form (11331) near the boundary, once we perform the change of coor- 
dinates^! 

2 _y 2 {l - a 2 cos 2 6-b 2 sin 2 0) + a 2 sin 2 6 + b 2 cos 2 - a 2 6 2 
2 sin 2 0= (?/2 + fl2)sin2 ^ (141) 



r sin 9 

r 2 cos 2 



2 :2n (y 2 + b 2 ) cos 2 9 



Retaining corrections only to order 0(l/r 4 ) relative to the leading order metric (I133p . the 



32 We thank Sangmin Lee for pointing out a typo in a previous version of this paper. 
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metric in our new coordinates become; 

ds 2 = - (1 + r 2 )dt 2 + + r 2 (d6 2 + cos 2 0d^ 2 + sin 2 0d0 2 ) + -^dr 2 + 2(m + f g) dt 2 

1 + r 2 r b Ag r 2 A^ 

2(2am + bq) sin 2 9 , 2(26m + ao) cos 2 , (2ma 2 + 26g) sin 4 6 1 , , 9 

-dtd0 — -5 dtdip 1 - ] -' 



r 2 Ag r 2 Ag r 2 A^ 

(2m6 2 + 2ag) cos 4 9,,, 2(2abm + a 2 a + 6 2 g) sin 2 cos 2 , , , , 
-# + 9A , d^d0 , 



r 



2 A 3 Y r 2 A 3 



where A# = 7 2 = 1 — a 2 sin 2 9 — b 2 cos 2 0. with 7 as in § §3.2 



(142) 



The unit normal vector to constant r slices is given by n r = -7= (with all other 
components zero). As our metric contains no terms that mix r with other coordinates at 
leading order , this is the same as rf = — 1 =. 

The Christoffel symbols (that are relevant for the calculation of the stress-tensor) as 
calculated from this metric up to the first subleading term (i.e. up to 0{\) terms) m r 
are given by 

t r f 4(m + abq) \ , _ 2{2am + bq) 



r 2 + lV r 4 (l-a 2 sin 2 0-6 2 cos 2 0) 3 y * r r 5 (l - a 2 sin 2 - 6 2 cos 2 0) 3 
^ _ 2(2bm + aq) , 2 (2am + bq) sin 2 



tr r 5 (l — a 2 sin 2 — 62 cos 2 0) 3 " r 5 (l — a 2 sin 2 — b 2 cos 2 

1 ^ 4(a 2 m + abq) sin 2 \ ^ t 2(2bm + ag) cos 2 



r 



r 4 (l — a 2 sin 2 9 — b 2 cos 2 0) 3 7 ^ r r 5 (l — a 2 sin 2 9 — b 2 cos 2 i 



r V> = 1 L _ 4(6 2 m + a&g) cos 2 \ ^ 1 



r 



r\ r 4 (l - a 2 sin 2 - 6 2 cos 2 0) V 9r r 

(2abm + a 2 q + b 2 q) cos 2 ^ [2abm + a 2 g + 6 2 g) sin 2 



^ r r 5 (l - a 2 sin 2 - 6 2 cos 2 0) 3 0r r 5 (l - a 2 sin 2 - b 2 cos 2 0) 3 ' 

(143) 



33 The leading behaviour at large r of the pure AdS metric given in terms of and rda where a 
represents the coordinates 9, (f>, ip, t. In the expression below, we have retained all coefficient terms of these 
'forms' that are at most of order \. 
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The extrinsic curvature, 0^, is given by 



i 



6* = - ( 1 + 1 ) 5 ( 1 + ^(3 + a 2 sin 2 (9 + b 2 cos 2 ^ 4 ^ 7 



R 11/9 / 1 m 7 /o 2 • 2 /j 72 2 fl , 1 n\ , 4afo 97 sin 
0£ = — a/ 1 + 1/r 2 1 — (3a sin 9 — b cos + 1) + 



^4 v / 



J = - v^TlF (l - ^(36 2 cos 2 (9 



> . ■■> „ , \ Aabqj 6 cos 2 9 
a 2 sm 2 + 1) ] + ^ 



, _ 2(2am + 6g)7 6 sin 2 2(2am + 6g)7 6 (144) 

<y<h — ~* <j+ 



r 4 * j.A 

; _ 2(2bm + ag)7 6 cos 2 9 ^ _ 2(2bm + aq)rf 



^ _ 2(2abm + b 2 q + a 2 q)rf sin 2 9 , _ 2(2bam + b 2 q + a 2 q)~{ 6 cos 2 9 







where 7 2 = - — , . $1 ,, — 5-5 ■ Therefore 

' 1-a 2 sm 2 8-b 2 cos 2 8 



= 0^ = 4 + ' 



It is easily verified that 0^ is traceless when the r dependent divergent terms are cancelled 
by the counter terms at the limit r going to infinity. After cancelling the divergent terms 
and then normalising it according to (11341) the stress tensor is given by 

II" = ™ (/(3 + a 2 sin 2 9 + b 2 cos 2 9) - ^ 7 6 ^ 
87rG 5 \ m J 

m (.A(a-.2 in 4ab 1 6 



7W - 1) " —f 



811G5 \ m 
- _ m f 6 f 3a 2 sin 2 9 - b 2 cos 2 9 + 1 \ Aabq < 

= v v ^0 ) ~ ^r 7 

m f 4 / 2 2 1 \ 4a6g 
7 M7 a + — q— 7 



87rG 5 \ \ sin 9 J to 

n </></> _ m /^e Z' 3o 2 cos 2 - a 2 sin 2 + 1 \ _ 4a6g ( 



87rG5 \ \ COS 2 /TO 

TO / 4 / 2 2 , 1 \ 4a& <? f 

7 M7 a + — T7T 7 



87rG5 \ \ COS 2 / TO 

n* = n* = ( (a - 2 *<) 7 6 n* = n* = ( (& - ^) 7 6 

= n" = (±L) (ab - 2 J^±3l) 7 e n - = -^ 7 4 - 



(145) 
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As we have explained in § §6.4[ q/m ~ 1/V+, and so all terms proportional to q in the equa- 
tion above are subdominant compared to terms proportional to m in the fluid mechanical 
limit r + — > oo. Dropping all q dependent terms, we find (11451) matches perfectly with the 
stress tensor as derived in (1301) and ( 1341) upon identifying (<px, 2 ) = (0, VO' (^1, ^2) = (a, &) 
and using (|9"4"1) . 

D Notation 

We work in the ( — h ++) signature, /x, z/ denote space-time indices, i,j = l...c label 
the c different R-charges and a, b = 1 . . . n label the n different angular momenta. The 
dimensions of the AdS space is denoted by D whereas the spacetime dimensions of its 
boundary is denoted by d — D — 1. In this paper we consider fluids on S D ~ 2 x R or 
equivalently S^ 1 x IR. Here we present some relations which are useful in converting 
between n, D and d: 

D = d+l = 2n + 2-(D mod 2) 
d = D-l=2n + (d mod 2) 



~d - r 






2 




2 



where [x] represents the integer part of a real number x. 

A summary of the variables used in this paper appears in table HJ 
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Definition 


Symbol 


Definition 


D 


Dimension of bulk 


d 


D — 1, Dimension of boundary 




iMewton Constant m Ado£> 


n 


[a/zj, no. oi commuting 


C 


no. of commuting R-cnarges 




angular momenta 


-RAdS 


Ada radius (taken to be unity) 
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Volume of S d , 




Horizon radius 


^mfp 


Mean free path, 


S 
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P 
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s 


Fluid entropy 


S 
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T 


Fluid temperature 




Pressure 


n { 


Fluid R-charge 




Proper R-charge density 




Fluid chemical potential 


r 


Thermodynamic potential (1I5T) 


$ 


Thermodynamic potential (jTTT) 


h 


V/T d , see flU 


hi 


dh/dvi 


v% 






Stress tensor 


■ J s 


Entropy current 


Jf 


R-charge current 




dx^/dr = 7(1, u), fluid velocity 




Angular velocities 


7 


^ 1/0 

(1 - v 2 y 1/2 


v 2 






Projection tensor, + u^u u 




see QZJ 


C,v 


Bulk, shear viscosity 




Heat flux, see (|HJ) 


K 


Thermal conductivity 


it 


Diffusion current, see (JTj) 




Diffusion coefficients 


E,H 


Total energy (125ft 


S 


Total entropy (126ft 


L a 


Angular momenta (1251) 


Ri 


Total R-charges (12B|1 




Angular velocities (139|) 


T 


Overall temperature (139|) 


Ci 


Overall chemical potentials (T39l) 




Partition function (|38f) 




Extrinsic curvature 


if" 


Boundary stress tensor 


n n 


Unit normal to boundary 




Induced metric of boundary 


T 


Integration constant (1281) 




Thermodynamic parameters (1661) 


m, q, s u Hi 


Black hole parameters 




1/(1 + «<) CT) 



Table 1: Summary of variables used. Numbers in parentheses refer to the equation where 
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